1.3 — Escolha do par de conectivas conjugadas, 
aplicável ao grafo dado 


a) Grafos totalmente ordenados; 

b) Grafos parcialmente ordenados; 

c) Verificação das diversas propriedades; 
d) Resumo e conclusão de 1.3. 


2.º Capítulo 


OPERAÇÕES SOBRE CONJUNTOS VAGOS 
DEFINIDOS NUM RETICULADO IN- 
TEIRO 


2.1 — Reunião/Intersecção de dois ou mais con- 
juntos vagos. 

2.2 — Composição de Relações Vagas. 

2.3 — Complementar dum Conjunto Vago. 


3.” Capítulo 


DESCRIÇÃO DOS DIVERSOS BLOCOS DO 
PROGRAMA 


3.1 — Fluxogramas;:. 

3.2 — Introdução do algoritmo para teste da exis- 
tência de circulações. 

3.3 — Anexos 


4.º Capítulo 


EXEMPLOS DE APLICAÇÃO DO PRO- 
GRAMA 


f.º Capítulo 


CARACTERIZAÇÃO DO CONJUNTO ONDE 
UMA RELAÇÃO DE ORDEM FOI IM 
POSTA HI 


1.1 — Introdução 


Seja dado um conjunto suporte O = (1,2,3,...M) 
com Card (0) = M. 

Seja IG a matriz quadrada mXm que descreve o 
grafo da estrutura ordenada (OQ, >), com elementos 
inteiros, obtidos a partir da seguinte aplicação: 


E or ÃO eg (0,1) 


O <= não existe o arco 


bb mosto E 
o” (iu. 0? =— : 
V (us, to) E O$, 8 (10, 102)  |I<=existe O arco 

WU ——— > w 


Define-se ainda uma matriz Xx == | Xi; | tal que 


mxm 


x=IG+I onde I é a matriz identidade, e ainda o 
conjunto G definido por: 


G = (um, 1) Q 4 (191, 102) E 1) 


1.2 — Verificação: das propriedades do grafo 
a) Existência e unicidade de ínfimo e supremo 


Ínfimo definido pela condição 


4 o 
É keo i=1 et, |) =0 


Supremo definido pela condição 


l 


E 
keo 


M 
» g(6,k)=0 
|=1 


b) Existência de circulação 
b.l — Circulação do tipo (1,1) (Reflexiva) 


Deve-se verificar g(1, 1) = O, Vieo 


LO 


Fig. 1 


b.2 — Circulação do tipo (1, j), (5,1) 
(Simétrica) 


Deve-se verificar: 
1] Ag (, J) = 


=1=>g(,D)=0V.. 
g(), 1) G Deo? 
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a Fig. 2 


A maneira prática de verificar esta condição é 
utilizar a sua negação dada por: 


Eh pe? "Lj, D=DAgO,DALO0 


b.3 — Circulação do tipo cadeia 
UR SR 
(Anti-transitiva) 


Deve-se verificar: 


g(i,) = 1 Ag G,k) = 
=| A...Aag(m,n=1=> N, J 


=> g(n, D)= | 
a 


Fig. 3 


0,7 (, j,...m, neo? 


onde p designa o número de pontos. 


Tal como na verificação de b.2, não se testou 
directamente esta condição. 

A sua verificação para todos os pontos do grafo 
exigiria a utilização dum espaço de memória dema- 
siado vasto, de modo a cobrir todas as combinações 
possíveis. 

O método seguido consiste em, partindo do ínfimo, 
verificar, subindo ao longo do grafo, da existência de 
pontos donde partam «caminhos» cujos pontos ter- 
minais não sejam o ínfimo. Esta operação faz-se até 
atingir o supremo. Veja-se o 3.º capítulo para análise 
do fluxograma respectivo. 


c) Existência de ramificação 
Para um dado nó K: 


Verifica-se a existência de arcos que partam de K 
através da condição: 


eep=)j=LosM 


definindo-se C, = Card((jeO : g(k,)) = 1)). 


Há ramificações se C, > 2 
Há unicidade de caminhos se C,= 1 
K é o ínfimo do grafo se €, = 0. 


O grafo diz-se totalmente ordenado quando se 
verificar a condição: 


E 4d 


k Ykeq 

Esta verificação faz-se partindo do supremo, até 
atingir o ínfimo, e à medida que se desce no grafo 
são construídos dois vectores que, no caso de grafos 
totalmente ordenados, são utilizados no cálculo do 
Máx/Min dum conjunto de pontos. 

São eles AUX (1) que contém os elementos do 
grafo ordenados segundo a relação de ordem im- 
posta, e INVAUX (i) que representa a aplicação 
inversa da que é definida por AUX (i). 


Exemplo: 
Supremo == V 
== > 


AUX (1) = V 
INVAUX (V) = 1 


d) Resumo e conclusão de 1.2 


Podem-se descrever as sucessivas verificações de 
propriedades do seguinte modo: 


Não existe ou não há unici- 
dade de infimo ou supremo. 
Não se trata dum reticulado. 


Há circulações. 
Não se trata dum reticulado. 


Não há ramificações. 
Grafo totalmente ordenado. 


Há ramificações. 


Grafo parcialmente ordenado. Fig. 4 
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1.3 — Escolha do par de conectivas conjugadas, 
aplicável ao grafo dado 


Vejamos agora quais as propriedades de que os 
pares de conectivas, 


(Máx, Min) ==(y, t) e 

(Sup, Inf.) ==(V, A), gozam para os seguin- 
tes tipos de grafos, que são os mais comuns 
e que esgotam as aplicações correntes. 


a) Grafos totalmente ordenados 


Neste tipo de grafos as conectivas (y, À) 
confundem-se com (V, A) pois quaisquer dois pon- 
tos são comparáveis entre si. 

No entanto o cálculo do Máx/Min torna-se mais 
rápido fazendo uso dos vectores AUX (i) e 
INVAUX (1) já referidos. 

Fazendo corresponder a cada tipo de grafo um 
valor diferente duma variável do programa (RETIC), 
pode-se identificar o tipo de conectivas a utilizar com 
o fim de se obter maior rapidez de cálculo. Ao 
mesmo tempo podem-se inferir certas propriedades 
comuns a todos os reticulados dum dado tipo. 

Dentro dos grafos totalmente ordenados podemos 
ter: 


a) Tipo RETIC 1 (M=2) 
RETIC = 1 
( M=2 
((0,1), 4, 4) é o Reticulado de 
0 Bool já estudado na Bibliografia 


referida. 
Fig. 5 


Os conjuntos vagos definidos sobre este reti- 
culado degeneram em conjuntos usuais. 
a;) Tipo RETIC 2 (M > 2) 
RETIC = 2 
M >2 
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É imediato por simples observação 
verificar as seguintes propriedades: 


(a) — Existência e unicidade de Ín- 
fimo e Supremo do grafo 

(b) — Fecho 
(c) — Unicidade no resultado 
(d) — Idempotência 
(e) — Comutatividade 
(1) — Associatividade 
5 (g) — Transitividade 

(h) — Absorvência 

(1) — Distributividade de y em re- 
6 lação a fh e de À em rela- 
lação a y 


Fig. 6 
(7) — Não complementaridade para todos os 


pontos à excepção do ínfimo e do su- 
premo que são sempre complementares. 


b) Grafos parcialmente ordenados 


Neste tipo de grafos já não se verifica a proprie- 
dade de fecho para as conectivas (y, 4) pois há ele- 
mentos que não são comparáveis entre si. 

Também dentro deste tipo podemos ter: 


b;)) Tipo RETIC 3(M=4) 


À RETIC = 3 
M=4 


Note-se que no caso do grafo 
da fig. 7 não existe y (2,3) 

b mas podemos ter V (2,3)= 
=|,e A (2,3)=4. 


Fig. 7 
São verif.cadas as seguintes propriedades: 


(a) — Existência e unicidade de infimo e su- 
premo 

(b) — Fecho 

(c) — Unicidade no resultado 

(d) — Idempotência 

(e) — Comutatividade 
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(f) — Associatividade 

(g) — Transitividade 

(h) — Absorvência 

(1) — Distributividade de XY emrelaçãoa A e 
de A emrelaçãoa V. 

(j) -— Complementaridade de todos os pontos: 


VA4,D)=1,A(1,9)=4 


b;) Tipo RETIC 4 (M>4) 


RETIC = 4 
M>4 
j | 
7 3 
2 5 
4 
3 
6 
7 
a) b) 
Fig. 8 


Qualquer destes grafos pode ser a reunião de 
diversos ramos totalmente ordenados com 3 ou mais 
elementos e com um mesmo supremo e ínfimo. 


b2.] — Reservando para o caso do Reticulado de 
Zadeh o termo «involução» que designa a 
diferença de um dado elemento para o 
supremo, define-se pseudo-complementar 
de x, como o (ou os) elementos (s) y tal 
que se tenha: 


V(gxy=VAA(,W=O 


Onde V designa o supremo e O o ínfimo do reti- 
culado. 

Não se exige portanto a unicidade. 

A multiplicidade de pontos pseudo-complementa- 
res dum dado ponto K, é função do número de 
ramos do grafo que não incluem o ponto K, e do 


número de pontos existentes em cada um desses 
ramos. 

Qualquer dos reticulados da figura 8 verifica a 
condição: 


Vecn d veo . V (x, y)=V A A (x, y) =) 


Um reticulado que goze desta propriedade, diz-se 
pseudo-complementado. 


b2.2 — A propriedade da distributividade poderá 
ser verificada ou não (por exemplo o grafo 
da fig. 8a) não verifica enquanto que o 
da fig. 8b) já verifica). 
O programa de computador verifica 
portanto se é satisfeita a condição: 


“itjeno SB AG] =A[V GV GM) 


e também 


“atear A VGI]=V [A 6,9, 46,0] 


Dado que se verifica a comutatividade, não é 
necessário ver também a distributividade à direita. 

Os grafos do Tipo RETIC 4 (fig. 8) munidos 
das conectivas Sup. e Inf., constituem reticulados 
que gozam as seguintes propriedades: 


(a) — Existência e unicidade de ínfimo e supremo 

(b) — Fecho 

(c) — Unicidade no resultado 

(d) — Idempotência 

(e) — Comutatividade 

(f) — Associatividade 

(g) — Transitividade 

(h) — Absorvência 

(1) — Pseudo-complementaridade (a estudar futu- 
ramente) 


b;)) — Tipo RETIC 5 (M>4) 


Neste tipo de grafo já há pontos que não têm 
complementar nem pseudo-complementar. 

(Exemplo: 3,5, 7). 

A distributividade é verificada pelo computador 
Também se determinam os pontos que têm comple- 
mentar. 
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RETIC =5 
M >4 


Fig. 9 


As propriedades verificadas são as seguintes: 


(a) — Existência e unicidade de ínfimo e supremo 
(b) — Fecho 

(c) — Unicidade no resultado 

(d) — Idempotência 

(e) — Comutatividade 

(f) — Associatividade 

(g) — Transitividade 

(h) — Absorvência 


by) Tipo RETIC 6 (M>4) 


Neste caso deixa de haver unicidade no resultado. 

Por exemplo W (7,8) tanto pode ser 4 como 5; 
ou A(5,6) tanto pode ser 8 como 9. 

Já existem subrotinas para tratar este tipo de 
grafos que, por serem susceptíveis de ainda sofrer 
alguns melhoramentos, apenas serão apresentadas 
futuramente. 


RETIC = 6 
M>4 


Fig. 10 


c) — Verificação das diversas propriedades 

O programa está protegido contra o tra- 
tamento de casos degenerados, sendo inter- 
rompido quando fôr caso disso, e imprimida 
a respectiva mensagem. 

Estas verificações constituem uma suces- 
são e cada passo é executado se os passos 
que o antecedem forem satisfeitos. 


(a) — Existência e unicidade de ínfimo e supremo 
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A não verificação desta propriedade im- 
plica a interrupção do programa, sendo 
escrito: 

«NÃO EXISTE O ÍNFIMO DO 
GRAFO» ou 

«NÃO HÁ UNICIDADE DO ÍNFIMO 
DO GRAFO», havendo mensagens idênticas 
em relação ao supremo. 


(b) — Fecho 
A não verificação desta propriedade im- 
plica a não verificação de (a), e o programa 
é interrompido. (ver anexo 1). 
(c) — Unicidade no resultado 
Verifica-se quando da chamada das subro- 
tinas Sup. e Inf., para a construção das 
matrizes MAI e MEN. Caso não haja unici- 
dade no resultado o programa é interrom- 
pido, e escrita a mensagem: 
«NÃO SE TRATA DUM RETI- 
CULADO POIS NÃO HÁ UNICIDADE 
DO RESULTADO COM OS OPERADO- 
RES SUP. e INF.» 


(d) — Idempotência 
A não verificação desta propriedade im- 
plica a existência de circulações de tipo 
reflexivo (1, 1). 
O programa é interrompido e é emitida 
a mensagem: 
«EXISTE UMA CIRCULAÇÃO NO 
PONTO K» (ver anexo 2). 


(e), (f) — Comutatividade e associatividade: 

A sua não verificação implica a exis- 
tência de circulações de tipo simétrico ou 
anti-transitivo, dado que neste passo já foi 
verificada a unicidade no resultado. 

Nesse caso é escrita a mensagem: 

«EXISTÊNCIA DUMA CIRCULAÇÃO 
NOS PONTOS I E J», ou então: 

«EXISTÊNCIA DUMA CIRCULAÇÃO 
QUE ENVOLVE MAIS DO QUE DOIS 
ELEMENTOS E CONTÉM O PONTO ko». 
(ver anexo 3). 


(g) — Transitividade 
Verifica-se sempre. (ver anexo 4). 


(h), (1) —- Absorvência e distribuidade 
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São testadas em conjunto, dado que à Pode definir-se um subconjunto de U através dum 
absorvência é um caso particular da distri ' caracterizante que será um vector com dimensão 
butividade, em que dois dos elementos rela- igual a Card(U). Cada componente representará o 
cionados são iguais. grau de pertença de cada elemento ao subconjunto 

de U assim definido. 
()) — Complementaridade 


Determinam-se sempre os pares de pontos Exemplo: 
complementares. 
Mesmo que o reticulado não seja com- U =(a,b,c,d,e,f) 
plementado, existe pelo menos um par de 
pontos complementares, que são o Supremo  Caracterizante de U = r,= (bode f V 


e o Ínfimo, podendo eventualmente existir 


mais pontos complementares, como no caso | | 
Caracterizante de q =T = (sbcdef) | 
p 


da fig. 11 em que 4 e 6 também são com- 000000 2 
plementares. 
Sendo A = (a/,, b/;,c/,, d/0,€/2, 1/1) 
Õ 
Caracterizante de AT es pa a] Fig. 12 


ou simplesmente T, = (v,1,1,0,2,1) 


Fig. ao ? oi o 
pd quando não há ambiguidade em relação à ordem 
pela qual se colocam os vários elementos. 
d) — Resumo e conclusão de 1.3 
S (sim) N (não) E Ullbece wii E | — Reunião /Intersecção de dois ou mais con- 
juntos vagos. 
Tipo ad a l 2 | 3 4 5 Sejam dados os conjuntos A,, A;, As,----, A, e à 
( ) | É | ; sua reunião e intersecção respectivamente: 
n 
Propriedades Ss S |:8 s 5 R=UA=AUAUA;U----UA, 
a) a h) i=1 
Distributividad Elm |? E 
O PISCA. N=NA=ANANAN----NA, 
i=1 
| Complementaridade S N S ? E 
| a com: 
Pseudo complementari 
dade S N S » ? = (&1, Bi, Y1,= ===, M) 
1 : | 
O Tipo RETIC 6 não se trata dum reticulado e por isso Es (az, Bay Y2, = ---, M2) 
fica fora deste quadro. : 
ae Capítulo A == (£,, Ba, à E ale Th) 
OPERAÇÕES SOBRE CONJUNTOS VAGOS O caracterizante do conjunto reunião será dado 
DEFINIDOS NUM RETICULADO INTEIRO UH] por: 
Sejam dados um universo U e um reticulado inteiro Potes DS Vocee=[1) 


(2, V, A). Seja dada uma aplicação T:U is É, onde 
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qa = SUP (m,q,,03,----- » La) 
8” = SUP (8:, B2,B3,----- » Ba) 
nº = SUP (m, É RS nedadiand , No) 


e do mesmo modo: 


ae (a, B, à Pa dia al ; M ) 
onde 
q = INF (a, Ga, Gan , E) 


B* e INF (8,, Ba, Ba, ac , B,) 


nº = INF (m, 112 03," Th) 


As expressões anteriores sugerem uma arrumação dos 
n caracterizantes em matrizes de dimensão nXcard 
(U), onde cada linha representa o caracterizante dum 
dado conjunto e cada coluna diz respeito a um dado 
elemento do universo. 

Para calcular o caracterizante do conjunto Reu- 
nião /Intersecção bastará aplicar o operador Sup. /Inf. 
sucessivamente a todas as colunas da matriz assim 
constituída. 


2.2 — Composição de relações vagas. 
Sejam dados: 


— um conjunto X com card (X) = 4 
— um conjunto Y com card (Y) = m 
— um conjunto Z com card (Z) = n 
— um reticulado (0, V, A) 


Define-se uma relação vaga J, como uma aplica- 
ção do produto cartesiano XX Y no reticulado consi- 
derado, e pode ser representada por uma matriz 
txXm (com elementos pertencentes ao conjunto OQ). 

Sendo dada uma outra relação vaga Jg; definida 
em YXZ como uma matriz mXn, a composição 
Jaop==Jao Jg será uma matriz tXn obtida através 
dum produto matricial, em que A está pelo produto 
e V está pela soma. 

Por exemplo Jaop (Xs, Z2) será o 2.º elemento da 
3.º linha da matriz da relação composta, obtido pelo 
produto da 3.º linha da matriz de J, com a 2.º coluna 


da matriz de J,, ou seja, usando o par de conectivas 


(V, A). 


Jos (X3, 2) = VÍA [A(X3 yo), Jo (Yi, Z2)], 
, A [Ja (xs, y2), Ja (Ya, Z2)] ' 
) A [Ja (X3, Y3), Jp (ya, Z2)] E) 


» 


, A [Ja (X3, Ym), Ja (Y mo Z2)| ) 


O resultado com o par (A, V) obtém-se de forma 
idêntica por simples permutação dos símbolos. 


2.3 — Complementar dum conjunto vago. 


Sendo dado um conjunto vago qualquer 
À, com 


Pa — (x, 8, “Ea ô, Et = dies + n) 


existe complementar de A se todos os elementos do 
reticulado, que figuram no caracterizante, tiverem 
complementar. 

Pode portanto existir complementar dum conjunto 
vago definido num reticulado não complementado. 
Basta para tanto que no caracterizante respectivo 
apenas figurem elementos para os quais existe com- 
plementar. 

Considerando ainda o reticulado da fig. 12 e o 
conjunto A definido por 

R=(V,l, 1,0,2,1) o complementar de A será: 

TE =0, 2,4 Vo do): 


3.º Capítulo 


DESCRIÇÃO DOS DIVERSOS BLOCOS DO 
PROGRAMA 


3.1 — Fluxogramas 
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a) Programa Principal 
Chama as subrotinas GRAFO; UNIRES; 
DISTRB; COMPLM; REUNI; INTERS; 
CMPRSI; CMPRIS; CMPLMT 


INTCIO 


LEITUAA DO lelpicA pod 
FARA cOMGIDA PO 
PAD GA AMA [ICORA) > 


ILORA vp = 


FIM 


[Erico] 


" LEITURA DO NUMERO DE 


fONTOS do bRAFO E pô mbi- 
CápOs FARA ESCAITA DE HATLILES 


LEITURA PA MATALE 
LEPLESENTAT IVA mé GEAFO 


ESCAITA DO CABEÇALHO É DA 
HATRIZ AE PAESEMTRTIVA RO LAP 


FILTRO DE PLOPAIE- 
| bapra po tMAFO DO 


MLRES 


A ILADADE JO AMSSLTADO 
cos do GPRLAROLES | 
sur E NF (FILTROS [| 


Dis TRIRUT NIDA BE 


Dãs dPpLADOLE 
ZuP E 1d F FTA E 


SOMAM O 


CONPLEMEA TRA DA DE 


ESLRITADAS HATAIEES 
AMIDA JRESOR 2 MALIHEM 


LEITULA sã I4EF FALA 
ESLOLMá SA SEQUEMC A 
2 Pça 


CEF = 


LEUANÃO DE pós DU 
MS LDA JATO ir 


COMPOSIÇÃO DE aELA 
CÕES VAÇÃOS (SUP/ IMF) 


e EMPRIS 1 
COMPOSIÇÃO DE AguA- 
ÇÕES VARA S(uE | SUP) 


[ CMPLMT | 
COMPLEMENTA bum 
CONJUNTO VAGO 
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b) Subrotina GRAFO (M, IG, INFIMO, SUPRE, 


c) Subrotina EXUNIS (M, IG, INFIMO, SU- 
PRE, *) 
Verifica a existência e unicidade de ínfimo 
e supremo do grafo 


IND, NC, *) 

Serve de filtro ao grafo para evitar tratamento 
de casos degenerados. Chama sucessivamente 
EXUNIS; EXCIRC; RETLIN 


[ Etynis 
TESTE Da ANS MEME 

E uNILIDADE D. 
FIM E SUPREMO [| 


TESTE DA ELISTEML A | 
| DE LIACULAÇÕES [y 


E RETLIN 
TESTE E TAATAMEMTD 
DE RETCULADOS EM Lisnhd 
[TOTAL EA TE dapéua DOS 


ESCAITA né HENSADÉM 
ÃO já WAJICI DADE 
DO inFIMO 


pi] 
[NEIAG=3] 


FSCAITA ME MENSAGEM 
“AO ElISTE O 
AMÉM 


RETURNA 


E 
ISUP pia 
“ EMAITA DE HENSAGERN 
Excel e 
DO SUPAEHO 


RETUANA 


[4 ESCRITA se MENSAGEM | 
NÃO EXISTE O 
ESCAITA DE MENSASEM SUPREMO 


EXSTENCIA E UN ÍCIDA DE 
RETURN 


DE PRO E BUPESHO 
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e) Subrotina EXRENM (M, LA, K, LAK) 
Extrai do vector LA, o elemento na posição 


dy Subrotina EXCIRC (M, 16, INFIMO, *) 
Verifica a existência de circulações. Chama 
as subrotinas EXRENM; ANSRCH 


costando à esquerda os elementos sucessores, 
e colocando um zero na última posição. 


|| een | 


(CAraI=O] 


| LApije INTIMO 
| LEtsi— INTIMO 


[ EXRENM 4 
ESTRAL INFIMO é 
AENUMEAS IE | 


EXTRA IAM) DE IA 
É dENLNEAA JÁ 


ESCRITA dE HENSABER | 
AUSÊNCIA 8 CIRCULAÇÕES 
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K (LA (K)== LAK) é renumera o vector en | 


f) Subrotina ANSRCH (M, AI, IA, IB, IC, IG, 


IAB, MA, MB, MC, MAB) 

Encontra os antecessores de AI que perten- 
cem a IC. 

Chama as sobritinas EXRENM, SCSRCH; 
ACTIAB. 


ACTUALIZA OS PAES 
PEATINLENTES À JAR 
| ni 


RETURN 


61 


g) Subrotina SCSRCH (M, BI, MB, MAB, IB, |h) Subrotina ACTIAB (MB, MAB, IB, IAB) 
Actualiza a matriz IAB extraindo todos os 
pares (x, y) tais que y € IB, Y (x,y) € TAB. 
Coloca x em IB, caso x tenha deixado de 
figurar na primeira coluna da matriz IAB ao 
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IG, TAB, IV) 


Encontra os sucessores de BI. Se todos os 
sucessores de BI pertencerem a IB, BI é colo- 
cado na última posição não nula de IB. Se 
houver um arco BI -» [2 tal que |, € 1B, en 
tão o par (Bl, 12) passa a constituir a última 


linha não nula da matriz IAB. 


[sesmen | 


IB(M6) «— BI 


ser feita a operação anterior. 


ABI 2) — O 
Aje- JAB(16,4) 
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1) Subrotina RETLIN (M, IG. INFIMO, SUPRE, 
IND, NC) 
Verifica se se trata dum reticulado total- 
mente ordenado (em linha). Sendo caso disso, 
ordena os elementos criando os vectores 
AUX (1) e INVAUX(I), e escreve as proprie- 
dades de que o reticulado goza. 


poem | 


Let 
1 suUPRE 
Ant] 
INvVAUN(T) 4 


Le—L+4 J-J+41 
ICo-ic+s4 
AuXiL])e= 9 


q 


NC 4 
NICOMP (1,4) e gurar 
MCONP(s,2) e INTIMO 


ESCRITA DE MENSAGEM 
acer EM LIMA 


EMAITA pr MENSAÇEMS ESCRITA é MENSAGEM 
AETICULADO DE BodL METICUCLADO JÃO 
(comPLEnENTADO con enEnNTtaDdo 


|» Subrotina UNIRES (M, INFIMO, SUPRE, *) 
Verifica se há unicidade no resultado com ns 
operadores SUP ec INF, e escreve as proprie- 
dades de que o reticulado goza, inferidas a 
partir do estudo feito. 
Chama a subrotina MAIMEN. 


ESCAITA p MENSAGEM 
FAO PLEDADES INFEAÍDAS 


| Subrotina MAIMEN (M, INFIMO, SU- |m) Subrotinas MAIORS/MENORS (M, A, NM, 


PRE)(3 Ka 
Determina o resultado da aplicação de SUP/ Constroem as matrizes MAIOR/MENOR de 
/INF a todos os pares de pontos do reticulado. valores lógicos com: 


Chama as subrotinas MAIORS; MENORS; 


SUP; INF. T<=12] 
MAIOR (1, J)= 
[E <=1*] 
T<=I<) 
e MENOR (I, J)= 
F<=14] 


poe | 


EEZE 
EEIE 


3)s=, FALSE. 


E 


PTE 


[etes || 


DETEAMINA ÇÃO 
DA MATE MALOR[I a) 


[RERSAS 
DETERMINAÇÃO 
DA MATRA MIESOL(L (3) 


3-VMIL) 
L-L-2 
TievMiL) 
Ie vm(L-A) 


<> 


Nota: 1 > ] significa; | não é maior ou igual a ], o que 
não quer dizer que | <' ], dado que | e |] podem 
não ser comparáveis. 
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n) Subrotinas FUNCTION SUP/INF (M, NP) o) Subrotina DISTRB (M, *) [2). [41 p) Subrotina COMPLM (M, NC, INFIMO, SU- 
fa), [4] Verificação das propriedades de absorvência e PRE, IND) 2). (4) 
Determinam o SUP/INF dum conjunto qual- distributividade. Verifica a propriedade de complementaridade 


quer de NP pontos do reticulado. para todos os pontos do rericulado. 


ZE 
2: 
14 
om 
É 


MH 
PI PTE Ch-i> 
IPA ga 


E , | a EST ] EM- [a] *ANFIHO | 
MR TSE NS gre JA | com 


Doo RETURN ) 
VBLIJ-FALSE. Spiova]> E 


Jo 1414 


Lt 


NCOMP(L A) Q | 
NCOMP(L,2) =D 


HAI ==HALÇT,]) 
MALE MATE) 
 MadR e MALI) 
HAN 1J-=MEN(T O) 
Md Lie MENS LR 
Hut MENS) 


[E] 


REN(LAGRT> 
= RI (MN Madi) 
PMEN(MAI,Haik)= 
= MAIT MNE) 


| RETURM 


ESCAITA Da MENSAGEM 
DETICULADO COnPLEMENTADO 


[FONTOS EC DAFLE HENTARE 5 


| RETURN 


ESCAITA DE MENSAGEM 
NãO ABSOMVENTE 


RE TURNA 


NMDSTa= 4 | 


E ea 
ESCRITA DE MENSAGEM ESCAITA DE MENSAGEM 
VO(I)= TRUE ABSOANENTE E ABSOAVENTE E 
já: Ê DISTRIBUTIVO “NÃO DIS TAIBUTIVO 
E | + E La] | + L+ 1 | 
- À - | viLjae- 1 
| (RETURN) RETURN 
p UM NC E SUPLE V 
I>M E (sp II 
E 


SO RETIC ho 


4 
SUP e-4[]) A 
EE: Sur 6 | 
Ú RE TURN RETURN 
NE = D 
INDe— 4 IND=— dé 
ESCAITA né MENSAGEM | | ESCRITA RE MENSAGEM 

BETeuLaDo COMPLEMEUTAnO) |RETICULADO mão COMPLEMENTADO 
PONTOS COMPLEMENTREES PONTOS COMPLEME NITARE 5 
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q) Subrotinas REUNI/INTERS (M, NE NUM- 
CON) 131. (4) 
Determinam o caracterizante do conjunto reu- 
nião intersecção de NUMCON subconjuntos 
dum universo com NE elementos. 
Chama a subrotina SUP/INF. 


LEITURA DO NUMERO DE 
IVBCON JUNTOS CUIA 
AEUNTÃO 3€ Quien CAM 


LEITURA DO NUMERO os 
ELEMENTOS DO UNIVERSO 


LEITURA pa MATRIZ 
sOs CARACTEMEANTES 
DOS CONTUNTOS 


ESCAITA DO CABEÇALHO 
E DOS CALACTERICANTES 
DOS CONJUNTO 


ViDeXYy(1,3) 
FESTUM 


INDAEU(I) e SUP(M NUM CON) 


a 


ESCAITA DO CAACTE RI 
ZANTE DO CONSUNTO 
REUNIÃO 


RETURN 


r) Subrotinas CMPRSI/CMPRIS (M, CARDI, 
CARDJ, CARDK) 3). [4 
Determinam a composição de relações vagas 
através dos pares de conectivas (LIM.SUP,, 
LIM.INF.) / (LIM.INF., LIM.SUP). 
Chama as subrotinas SUP; INF. 


LEITURA 50 CARDINAL 
DOS CONTUNTOS 1,3, 


LEITURA DAS NATRIRTS 
DAS AELAÇÕES À s E CUJA 
conrgsição É "AE TENDEDA 


VIA (1,9) 
vi(2)—XB(3,4) 


ESCAITA so CABEÇALHO 
E 


COMPOSTA sE A com & 
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s) Subrotina CMPLMT (IND, INDCON, IND- 
COM, NC) f3). [41 
Determina o complementar dum dado con- 
junto vago, caso exista. 


| CMPLMT | 


LEITURA vo NÚMERO sé 
ELEMENTOS Do UNIVERSO 


ESCRITA DE MENSAGEM 

AE TICULADO mão COMLEMENTADO 

NÃO EXISTÊNCIA se COMPLEMENT 
»O CONJUNTO DADO 


RETURN 


ESCRITA E MENSAGEM 
RE Ticui NÃo COMPLEMENTADO 
HAS EXISTÊNCIA DE COMPLEMENTAR 


ESCRITA DO CARACTERICANTE 
Do COMPLEMENTAR DO 
CONJUNTO DADO 


RETURN 
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3.2 — Introdução do algoritmo para teste da exis- B =(B;) AB= (AB) 1.7 sob 

tência de circulações. 

com Ai, B,, C, CP (0) 
a) Conceito de operadores ANT(x) e AB, co? 
SUC (x) 
Satisfazendo as condições seguintes: 
Designa-se por ANT (x) o conjunto dos anteces- 
sores de x dado pela condição: | G = 09 


C, = CW, N ANTE» (INFIMO) k= 1,2,...,n 
ANT(x) = (ye c |g(y, x) = 1) 


A, = (xeQ |xe ANT? (INFIMO) A 
e por SUC(x) o conjunto dos sucessores de x dado 


apa AxeCijk=1,2,..,n 
pela condição: 
B, = INFIMO 
SUC(x) = (yeo|g(x,y) = 1) 
B, = (xeon|(xe By) V 


Simboliza-se por ANT a aplicação sucessiva do 


V (x e A, A SUC (x) B,1) V 
operador ANT aos elementos sucessivamente gera- 


V((x, ad YV 
dos por ele a partir dum certo elemento «a € 0. (6x 1) 6 Als A pel ) 
AB, = 
ANT (x) = Ae O e 
AB, = ((x,9) € G| ((x,y) AB, A 
ANTO) (x9) = ANT [ANT(x0)] = ANT(A;) 


Ay éB)Vxec A Ay é By)) 
= A=(yeQ|ye ANT(x;), VxcAi) 


ANTO (x)= A;= (yeQ|y e ANT(x;)), Vx;ceA2) A condição de não existência de circulação é 
dada por: 


= G=4 => B=0 ou ainda por 
ANTO (x)= A,=(yeQ e ANT(xa), VXn € Ar 
(Xo = (y [Y ) k-1) Udo A má 


Poder-se-á condensar mais a nomenclatura fazendo: 
SUCO? (x) = x = ANT? (x) 
SUC€” (x) = ANTC? (x) Exemplo 1 


SUCO (x) = ANTCD (x) 
ou ANT” (x) = SUCO” (x) 


ANT(CO (x) = SUCCD (x) Fig. 13 


b) Introduzem-se as seguintes sucessões de con- 
juntos: 


A=(A) C=(C) 
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ES C, A, B, AB, 
9) | soa e Bê 
0 |M,23,4,5,6,7) (8) o 9 
(1,2,3,4,5,6) | (7,3) (8,7) (5,5); (3,6)) 
(2,5) (1,4,6) 18746) | GA) | 
9 (2,5) (8,7,4,6,5,3,2,1) $ 
Exemplo 2 
10 Fig. 14 
i Cc, B, 
0 2. | — | — 
| 1 (1,2,3,4,5,6,7,8,9) (10) 
(1,2,3,4,5,6) (7,8,9) (10,8,9) ((7,3)) 
1.3 (1,2,3) (4,5,6) (10,8,9,6) ((7,3); (4,7); (5, 7)) 
- (1) (10,8,9,6) ((7,3); (4,7); (5, 7); 
(2,4); (2,5); (3,5)) 
5 $ (10,8,9,6) (07,3); (4,7); 6,7); 


(2,4); (2,5); (3,5); 
(1,2); (1,3)) 


3.3 — Anexos 
Anexo 1 — Suponha que não existe V (a, b) 
O que pode ser traduzido por: 


(xe92]V(a, b)=x)=64=>(xe0] 
IV(am)=x) AIV(b,xy)=x])=gó 
Esta última condição é incompatível com a exis- 
tência de supremo, logo garante-se a propriedade de 
fecho, uma vez verificada a existência de supremo. 
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Anexo 2-— A condição de não existência de cir- 
culações de tipo reflexivo (1,i) é dada por g(i, = 
= 0, Vieo. Portanto IG é uma matriz de diagonal 
nula. 

Dado que se faz X=IG+I, obtém-se em X uma 
matriz de diagonal unitária, impondo-se deste modo 
a idempotência. 

Anexo 3 — Suponha-se que: 


SUP (a, SUP (b, c)) - SUP (SUP (a, b), c) 
Faça-se SUP (a, SUP (b, c)) = x (1) 
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e SUP (SUP (a,b), )=y (2) 


de(1): e de (2): 
x2a y>za 
| x>b Iy>b 
Xx2c Yac 


Vejam-se os casos possíveis: 

(x>y)=> x não pode ser o supremo pois y<x e y 
também é majorante de a,b,c. 

(y>x)=> y não pode ser o supremo pois x<y e y 
também é majorante de a, b,c. 

Portanto a única solução é termos x==y, ou seja 
SUP tem de ser associativa. 

Para a comutatividade o procedimento é idêntico. 


Nota: 

A unicidade do resultado garante que se xLy € 
y<[x então x ==. 

A não existência de circulações garante que se 
Xx>y então ypbx. 


Anexo 4 — Suponha-se que se tem para 3 pontos 
do reticulado 


SUP (a, b) = a A SUP (b,c) = bA SUP (a,c) <a 


Já vimos (associatividade) que: 

SUP [SUP (a, b), c] = SUP [a, SUP (b, c)] 
dado que SUP (a, b) = ae SUP (b, c) = b, vem: 
SUP (a, c) = SUP (a, b) 

SUP (a. c) = a 


Portanto, a hipótese: SUP (a,c)5=a é impossível 
desde que se verifique a associatividade. 


4.º Capítulo 


EXEMPLOS 


Fig. 15 
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4.1 — O «output» que se segue diz respeito ao reticulado da figura 14. 


OPERACÕES COM CONJUNTOS VAGOS (FUZZY SETS) 


*“1I=4ATRIZ REPRESENTATIVA DO GRAFO 


SIGNIFICA OUE EXISTE O CAMINHO DE 1 PARA J 


-- 
Pd 
had 
= 
Do uw 
So am 
[= 
jo 


te 
[ii] 
di 
— 
o 
e 
o 
o GQ 


un E pá po + 
D20 
Do 
20 0 
Doo 
Dos 
Dos- 
e 
=20E 
Doo 
ooo 


owspin 
odio da a 
om o] 
om-be 
Doo 
sao 
ocop: 
sooob I 
| 
sco 
o00| 
pa quit 


2= PROPRIEDADES DD RETICULADO 


AJEXISTENCIA E UNICIDADE DE 1 INFLMO —- PONTO 10 
EXISTENCIA E UNICIDADE DE 1 SUPREMO = PONTO 1 


MAIS DO QUE DOIS ELEMENTOS E CONTEM O PONTO 3 


MATRIZ IABITod) 
a 


= 
Lil ho UU é) 


- 


4.2 — Neste caso trata-se dum reticulado totalmente ordenado. 


UPERACÕOES COM CONJUNTOS vaGos (FuZzr SETS] 


i=MATRIZ REPOESENTATIVA DO GRAFO 


STENÍEICA. QUE EXISTE O CAMINHO DE | PARA J 


i TE TE RE RR E: 
10000000» 
2 DO DO DADO 
D O 1 O É uq 
0 1000000040 0 8 
SD 0 6 606 1LãHHãoOS ú 
SA OD. OD | 
Ce MEM TT MES a 
8sB92 0000000 06 
3000006004 6 À 
10. 0 090 10.6 60 DD 0 6 


êé= PROPRLEDADES DO RETICULADO 


AJEXISTENCIA E UNICIDADE DE | IMFIMO = PONTO & 
EXISTENCIA E UNICIDADE DE | SUPREMO - PONTO 7 
“AD EXISTE QUAL QUER rio DEE LEGUL ACAO 


E HA UNICIDADE NO ME SUL TADO 


RAT É DUM RETICULADO EM LINHA E JU 
OS ELEMENTOS ESTAO ORDENADOS DO SEGUINTE MODO 
[Li 


d 
5 


—-€£€A«zãõm eee == 


5 


= 


E ie jm a a ida 


DIIDEMPOTENTE 
FICO4MUTATIVO 
FIASSOCLATIVO 


e eme MEABSMEMENTE ll DE 


E PD ES a qui 


Pd pn ndo pa 
ES Le) ME 


BEsTe “990 O RETICUL ADO NAO E COMPLEMENTADO, 
2015 05 DuUTãos PONTOS NAD TEM COMPLEMENTAR 
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4.3 — Neste exemplo é tratado o reticulado da fig. 15. 
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OPERACOES COM CONJUNTOS VAGOS (FUZZY SETS) 


I-MATRI Z REPRESENTATIVA DO GRAFO 


SIGNIFICA QUE EXISTE O CAMINHO DE | PARA J 


SE A A E qa fio = 
DD0/jpoDOdjDD = 
CoBDocopDe= ly 
2090 00090) + Jd 
0000 C0DDdDrD É 
32200000 [=D iu 
00608600 =D0 ik& 
DdoDd= DSO | 
2000-DOIDO É 
00» = 0DD id 
md 
—- e es Sao =] E 
| 


2- PROPRIEDADES DO JETICULADO 


AJEXISTENCIA E UNICIDADE DE 1 INFIMO - PONTO 10 
EXISTENCIA E UNICIDADE DE 1 SUPREMO - PONTO 1 
REA 7 A! JER TIP 

RJjJAS OPERACDÕDES SUP E INF SAD FECHADAS 
E HA UNICIDADE MO RESULTADO 

CIJIDEMPOTENTE 

EJASSOCILATIVO 

FJTRANSITIVO 

GIABSORVENTE 


PROPRIEDADES ADICIONAIS 


=— o "»"»" MM "hm .——————oooeeeeeoeeeeeee ass 22 mm 


HINAO DISTSLauTivo 
PONTE! E a (aa LE [| T 
OS PONTOS 4 E & SAO COMPLEMENTARES 


= E MOC LA MPLEM A 
POTS 05 DUTOS PONTOS MAD TEM COMPLEMENTAR 


MATRIZ MALOGIIod) 


E E E F F E 
F F F F F F F F T T 
F F F r F F F F F T 


“ATi? ALHgBVI,dI] 


T [a r Fr F ” F r F F 
T T r F r F ) F r F 

US an SECA Song T ' T F JF Fr a E F ' E É o Des o 
sam E A T | - T F F F Fr F F 
T T T F T F [o r F F 
no T F T F F T F F F F 

AS T T T T T F T E Tr F = 

DUDE ren 7 T T F T E. [7 T E F a 
T T T k T T F F T F 

REA T T T T T T T T T T ra 
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AMATI17 MALIIod) 


MATRIZ “HEN[Lod) 


à- DPERACOES SOBRE CONJUNTOS VAGOS CUJOS CARACTERIZANTES 
SAO DEFINIDOS SOBRE O RETICULADO CONSIDERADO 


A-DE TERMINAL A A REUNTÇÃO DE 3 SUBCON JUNTOS 
DUM UNIVERSO COM 5 ELEMENTOS 


CADA COLUNA CORRESPONDE A UM DADO 


NOTA: CADA LINHA REPRESENTA O DC RR ra RA DE TON TVE 


TO DO UNIVERSO 


ELEMENTOS 1 2 E | & 5 
1 de Bo " 


CAJALTTERIZANTE DO CONJUNTO E: ( De Ze 6. do 5) 


O CARACTERIZANTE DO CONJUNTO REUNIÃO E DADO POR 


“H=-DETERMINAÇCÃO DA INTERSECCADN DE & SUBCONJUNTOS 
DUM UMIVERSO COM 5 ELEMENTOS 


&I A EL L: à DO E 
DOS NUAITS SE QUER DETERMINAR 


Lo = 


A INTERSECCAD. 


=SLEMENTUS É É So a Is 
CARACTERIZANTE LO COSJUNTO 1: [de Sa Be Told) 
CAJACTFERIZANTE DO CONJUNTO 2: [0d Ze 6 do 5) 
CARACTERIZANTE DO CUNJUNTO 3º [Te 3aLO, d+ 9) E 
CAJACTERIZANTE DO CONJUNTO q 


ia ia Mo CARATTERLIZAMTE. mM: JC is it Ti IMTE ESECCADn FF FNetal POB. ses 


CARACTERI Z. DU CONJ. INTERSECCAD: (1%, TelD. 7.10) 
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7a 


QMINACÃO DA COMPOSICAD DE DUAS RELACDES VAGAS 
ZANDO OS OPERADORES (LIM. SUP. É LIM: INF.) 


rim 
tem 


LAÇAO à E UAáA APLICAÇÃO DO PRODUTO CARTESTANO [ãJde NO RETICULADO CONSIDERADO. 
U MUMERO DE LINHAS = TGUAL AD CARDINAL DE 1 

U NUMERO D= COLUNAS E IGUAL AM CARDINAL DE J 

COMO O CARDINAL DE 1 E & .E O CARDINAL DE JE 4 MATRIZ DA RELACAO d 
VAL SER UMA MATAL? DE DIMENSADO CARDIZCARD J,OU SEJA UMA MATRIZ 4% 2a 


EE 


MATRIZ Da RELACAO à 


E N— Em EEE EE 


dy =X 

. l 3 ad 
LO a 

E e RE 
& T 


À RELACADO B E UMA APLICAÇÃO DO PRODUTO CARTESIANO déKo NO RETICULADO CONS LIDERADO. 
d NUMERO DE LINHAS E IGUAL AD CARDINAL DE 4 E 

PALM E LUNI UAL | CA [ 
COMO QU CARDINAL DE 4 E 2? ,E O CARDINAL DE K E 5 sã MATRIZ DA RELACAO B 
VAL SER UMA MATRIZ DE DIMENSÃO CARDJECARDE OU SEJAsUMA MATRIZ 2% Se 


à RELACAO COMPOSTA = UMA APLICACAD DO PRODUTO CARTESIANO [$4K+. NO RETICULADO CONSIDERADO. 
Jd NUMERO DE LINHAS E TGUAL AO CARDINAL DE 1 . 

U NUMER L | DINAL D K 

A MATA? DA RELACADO COMPOSTA VAT SER PORTANTO UMA MATRIZ DE DIMENSAD CARDISCARDE 

OU SEJA UMA MATRT] E bs Sa 


MATRIZ DA RELACAD DATIDA à PAR R Ê - 


DA COMPOSICAO DAS RELACOES à E B 


E? fs 2 LO & 


7 4 a 10 7 
Y-DE TERMINACÃO DA COVPOSICAO OF DUAS RELACDES VAGAS 
UTILIZANDO OS NPERADNAFS [LIMetNF. q LIM SUP. | 
À MELAÇÃO & E UMA ANLICAÇÃO DN PRODUTO CABTFS[aNO 144, NO SETFICULADO CONSIDERADO. 


do SUMETO Ci LINHAS = TGOJAL AU CARDTMHNAL DE T 

Jd HUMERO DE COLUNAS d LGUAL AM CARDTHAL DE 4 
4 4 Pre) E E E IRIZ.BA GELACAD À 

val vER UMA MATO? 96 DIMENSAMN CAGÓ9TSCARDO, OU SEJA qUMA MATAIE dv das 


MATALZ DA RELACAO A 


19 l E A 


a q da 


o Ra PS ea BA PE E Er Lo PAL 

JS NUMERO DE LINHAS E LúyAL AU CARDINAL DOE q 

à AUMESNTD 902 COLUNAS E TodJár 40 CARDIMAL DE 4 

Id O CASDImáL DE do e & 4 DO CAIDINAL DE &K E 5a vã MATRIZ DA RELACAD BR 
] L LT ' Fa FP [= MM ld E [e] ) | T : 


MATILr DA PELAÇÃO H 


a — ai sit POÇÃO 2 e essi eta 
e) 7 f 1 5 
Li a LO E” q 
A E [e 2 70 ê . 


A KELACAU COMPUSTA 2 UMA APLICACAO DO PRODUTO CASTES [ANO [is NO RETICULADO CONSIDERADO, 
do NMUMIRO DD: LINHAS 2 [QUAL AO CARDINAL DE 1 
E S NUMERO DE COLUMAS E IGUAL AO CAPDINAL SE k 
A 4ATAis VA RELA - U- 
DU SEJAcUMA MATRI? do Ga 


Ta val SER BOSTANTO UMA VATERLF DE GIMENSADO CARDIZCARDKE 


LO E E LO 
TS mão ? q = 
E Fê E E pj 


E -— 
La E EEE E ——— 
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E-DETEPMINAÇCAO DON COMPLEMENTAR DUM DADO CONJUNTO VAGO.CASO EXISTA 


CARACTERIZANTE DE 4 = (4,4. 3, Bo TelO. Do 


2. e Le De To 


3.10, 40, De. De de, 6 2.10. Be1lO, 6. 4. 8) 


Seo 


OD RETICULADO NAO E COMPLEMENTADO, 
“o 
PURTANTO NAO EXISTE COMPLEMENTAR DO CONJUNTO DADO 


À 


a Es E Em 


TF -DETECMINAÇAO DO COMPLEMENTAR DUM 


CARACTERNIZANTE DD à = (4 44 44 Se bo To 


O. 6.10+ 


DADO CONJUNTO VAGO,CASO EXISTA 


QGuelD. Le de GelOe de Oo De de. Te 1010, 6 4. 6. 6) 


E 


APISAR D? PECTICULADO NAM SER COMPLEMENTADO, 


FXISTE CIMPLEMENTAS DO CONJUNTO DADO POIS NO SEU CARACTERIZANTE 


APTNAS FIGURAM PONTOS 9UF TEM COMPLEMENTAR 


CARACTTRI ZANTE DN COMPLEMENTAR DO CONJUNTO DADO 


CARACTERIZANTE DO CoMP(Aa) = ( 6 
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RESUMO 


Apresentam-se as principais propriedades dos zeóli- 
tos. Os obtidos por síntese são muito mais utilizados do 
que os zeólitos naturais, pois devido ao grande avanço 
tecnológico das técnicas preparativas, conseguem-se ob- 
ter materiais com elevado grau de pureza e as proprieda- 
des mais convenientes para as aplicações industriais. 

Descrevem-se as suas principais aplicações industriais 
como adsorventes, em secagem, purificação de gases e 
separação de hidrocarbonetos, como permutadores ióni- 
cos, em recuperação de isótopos radioactivos, tratamento 
de esgotos e detergentes, e como catalisadores, em pro- 
cessos de refinação de petróleo e petroquímica. 


1. INTRODUÇÃO 


A origem do termo zeólito remonta ao século XVIII. 
Foi então constatado que determinados materiais quando 
calcinados a elevadas temperaturas pareciam, simulta- 
neamente, fundir e ferver com abundante libertação de 
vapor de água. Foram então designados por zeólitos, por 
composição das palavras gregas “'zein”” (ferver) e “li- 
thos” (pedra). 

Os zeólitos eram, então, considerados como simples 
curiosidades mineralógicas e, só na segunda metade do 
século XIX e princípios do século xx, começaram a sur- 
gir os primeiros trabalhos sobre as suas propriedades de 
permuta iónica e de adsorção. 

Em 1958, Eichorn verificou que os zeólitos naturais, 
como a chabazite e natrolite tinham a capacidade de per- 
mutar catiões sódio e cálcio de soluções diluídas. Wei- 
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gel e Steinhoff, em 1925, notaram que a chabazite ad- 
sorvia vapor de água e álcool etílico mas não a acetona 
e o benzeno, donde a designação dos zeólitos por pe- 
neiros moleculares, atribuída a McBain'», em 1932. 

A primeira aplicação laboratorial dos zeólitos desi- 
dratados para purificação de mistura de gases, deve-se 
a Barrer?), que em 1945 utilizou também a chabazite. 

No entanto, dado que os zeólitos naturais se mos- 
travam insuficientes para satisfazer os investigadores e 
industriais, em 1948 começaram a surgir os primeiros 
zeólitos obtidos por síntese, nomeadamente a morde- 
nite%) e os zeólitos A, X e Y 45.9, 

A partir de 1954, estes novos zeólitos sintéticos co- 
meçaram a ser utilizados como adsorventes e peneiros 
moleculares em processos industriais de separação e pu- 
rificação de gases, tendo actualmente inúmeras aplica- 
ções, com saliência para os processos de remoção de va- 
por de água e gás sulfídrico do gás natural e gás de 
petróleo liquefeito (L.P.G.). 

A primeira aplicação industrial dos zeólitos como 
permutadores iónicos, verificou-se na década de 60, na 
recuperação e concentração de estrôncio e césio radioac- 
tivos existentes em águas residuais alcalinas radioacti- 
vast), Actualmente, são inúmeras as aplicações em agri- 
cultura, nutrientes para animais, detergentes e tratamento 
de águas. 

A mais importante aplicação dos zeólitos observou- 
-se em catálise heterogénea, tendo o primeiro sucesso in- 
dustrial dos catalisadores zeolíticos ocorrido em 1962, 
no cracking catalítico. Os zeólitos vieram então, substi- 
tuir os catalisadores clássicos de silica-alumina amorfa, 
sendo a incidência económica dessa substituição parti- 
cularmente importante, nos processos de cracking diri- 
gidos para a produção de gasolina. 

Actualmente, os zeólitos devido às suas excepcionais 
propriedades de actividade, estabilidade e selectividade 
de forma'são muito utilizados nos processos catalíticos 
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da refinação de petróleo e petroquímica! ?, nomeada- 
mente no cracking, hidrocracking e produção de gaso- 
lina a partir do metanol. 


2. ESTRUTURA DOS ZEÓLITOS 


Os zeólitos são silicoaluminatos cristalinos, cuja es- 
trutura tridimensional apresenta uma porosidade regu- 
lar de dimensões comparáveis às das moléculas orgâni- 
cas, sendo as aberturas dos poros variáveis de 3 a 10Aº 
conforme o tipo de estrutura. 

Os elementos estruturais de base são os tetraedros 
AIO; | e | SiO, |, ligados entre si pelos quatro vérti- 
ces de oxigénio, originando assim, uma estrutura micro- 
porosa. As cargas negativas dos tetraedros | AIO; | são 
compensadas por catiões alcalinos, que podem ser subs- 
tituídos por outros catiões por permuta iónica. 

Conforme o arranjo dos tetraedros, assim se obtém 
uma grande variedade de estruturas cristalinas, conhe- 
cendo-se actualmente mais de trinta zeólitos naturais e 
uma centena de zeólitos sintéticos, embora só alguns 
como a mordenite, os zeólitos A, X, Y e os da série 
Z S M (Zeolite Synthetic Mobil) tenham interesse indus- 
trial. 


Fig. 1 Estrutura porosa do zeólito Y, com corte 


passando pelas diferentes cavidades 


Na figura 1, apresenta-se a estrutura do zeólito Y que 
contém 192 tetraedros, constituindo 16 prismas hexago- 
nais, 8 cavidade sodalite e 8 grandes cavidades. O corte 
apresentado na figura 1 permite melhor visualizar as li- 
gações entre as diferentes cavidades, podendo-se cons- 
tatar que as grandes cavidades têm 12,5 Aº de diâmetro, 
sendo as aberturas dos poros de 7,5 Aº. 

Ao contrário do zeólito Y, a estrutura da mordenite 
não contém cavidades (figura 2), apresentando grandes 
canais paralelos de abertura elíptica (6,7 x 7 Aº) ligados 
entre eles por pequenos canais paralelos (2,9 x 5,7 Aº). 
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Dado que a maior parte das moléculas orgânicas não 
podem penetrar nos pequenos canais, pode-se conside- 
rar a estrutura porosa da mordenite como monodimen- 
sional. 


(b) 


Fig. 2 Estrutura porosa da mordenite (a) 
e do zeólito ZSM-5 (b) 


Pelo contrário, a nova geração dos zeólitos sinteti- 
zados pela Mobil, como o ZSM-S (figura 2) apresenta 
dois tipos de canais de dimensões análogas, o que lhes 
confere uma estrutura tridimensional. 


3. PREPARAÇÃO DOS ZEÓLITOS 


Os zeólitos utilizados na indústria e na investigação 
são previamente submetidos a diversos tratamentos tér- 
micos e permutas iónicas que podem modificar as suas 
propriedades. 

Os zeólitos são obtidos por sintese hidrotérmica a 
partir dum gel aquoso alcalino de alumino-silicato, ocor- 
rendo a cristalização por um processo de nucleação con- 
trolada a temperaturas entre 80 e 300ºC. 

Uma grande variedade de materiais de silica e alu- 
mina, como silicato de sódio, gel de sílica, sulfato de alu- 
mínio ou alumina, podem ser utilizados para a obten- 
ção do gel, sendo o zeólito obtido, dependente dos 
materiais seleccionados. Também o tipo de catiões pre- 
sentes (sódio, potássio, etc.) tem influência na estrutura 
do zeólito. 

Após a síntese, os zeólitos têm os seus poros ocupa- 
dos por moléculas de água. Para utilização como adsor- 
ventes, é necessário eliminá-la por tratamento térmico 
(300-500ºC) sob corrente de ar seco. As cavidades zeo- 
líticas podem assim ser ocupadas por outras moléculas, 
o que confere aos zeólitos excepcionais propriedades ad- 
sorventes, pois a interface sólido-gás ou sólido-líquido 
encontra-se essencialmente no interior da rede cristalina 
porosa (volume poroso pode atingir valores de 50%). 

Os zeólitos, para utilização em catálise, devem pos- 
suir uma elevada estabilidade térmica. Por exemplo, os 
catalisadores zeolíticos de cracking, desactivados pela de- 
posição de coque, sofrem tratamento térmico de rege- 
neração sob ar (pressão de 2 a 3 atmosferas) a tempera- 
turas compreendidas entre 650 e 760º€C. 
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Para estabilizar a sua estrutura, cuja resistência de- 
pende essencialmente da natureza dos catiões de com- 
pensação e da razão sílica alumina, os zeólitos são sub- 
metidos a tratamentos térmicos a temperaturas elevadas, 
usualmente sob corrente de ar seco, que diferem con- 
forme o tipo de zeólito. A obtenção da forma ultra-está- 
vel do zeólito Y é conseguida por calcinação em presença 
de vapor de água, a temperaturas de 600-700º Co. 1D, 

As propriedades catalíticas dum zeólito dependem da 
natureza e teor dos catiões que podem ser introduzidos 
por permuta iónical2. 13), devido à excelente capacidade 
de permuta destes materiais. Podem-se permutar catiões 
metálicos que após calcinação e redução sob hidrogé- 
nio, originam centros metálicos que catalisam reacções 
de hidrogenação-desidrogenação. 

Os centros ácidos que catalisam reacções de cracking 
e isomerização podem ser obtidos por calcinação do zeó- 
lito na forma amónio, obtido por permuta catiónica com 
um sal de amónio. 


4. APLICAÇÕES INDUSTRIAIS 


As principais aplicações dos zeólitos são em adsor- 
ção, permuta iónica e catálise. 


4.1 Aplicação em adsorção 


Durante muito tempo, os adsorventes usados na in- 
dústria foram os gels de sílica ou alumina e carvão acti- 
vado. Devido às exigências da pureza dos produtos para 
a investigação e indústria, tornava-se imperioso desco- 
brir novos materiais, tendo os zeólitos um enorme su- 
cesso em adsorção devido às suas propriedades particu- 
lares que poderemos assim resumir: 


— selectividade geométrica, que permite aos zeóli- 
tos actuarem como peneiros moleculares 

— condensação capilar nos poros devido às suas pe- 
quenas dimensões, aumentando-se a concentração 
superficial do adsorvido 

— selectividade energética de adsorção, resultante da 
existência de cargas eléctricas na estrutura, mais 
ou menos separadas, que conferem uma grande 
afinidade para moléculas polares e polarizáveis. 


Os zeólitos mais usados em adsorção são sintéticos 
do tipo A e X, mais ricos em alumina e, consequente- 
mente, com maior concentração de catiões de compen- 
sação e maior selectividade energética de adsorção. 

As principais utilizações industriais são em secagem, 
purificação de gases e separação de hidrocarbonetos. 


4.1.1] Secagem de gases 


Os zeólitos apresentam em relação aos adsorventes 
clássicos (sílica-gel e alumina) uma capacidade de ad- 


sorção extremamente elevada para água e compostos po- 
lares com baixas concentrações na fase gasosa, como se 
pode verificar pela figura 3. 
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Fig. 3 Isotérmica de adsorção de água para três tipos 


de materiais 


Embora a capacidade de adsorção diminua para tem- 
peraturas mais elevadas, em zeólitos tipo A e X essa ca- 
pacidade é muito menos afectada, comparada com os 
adsorventes clássicos (figura 4), o que permite uma se- 
cagem eficiente para gamas de temperaturas entre 100 
e 150ºC€, conseguindo-se reduzir os teores em água a 
Il ppm para gases e 10ppm para líquidos. 
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Fig. 4 Isobáricas de adsorção de água para três tipos 
de materiais 


Com o zeólito 3 A, só as moléculas de água ( 34º) 
e amoníaco podem penetrar nos seus poros, conseguindo- 
-se, assim, uma adsorção selectiva da água. Isso explica 
a sua elevada utilização em secagem, embora não tenha 
uma grande estabilidade térmica. 
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Como algumas das principais aplicações industriais 
em secagem podemos referir as seguintes: 


— secagem de ar, hidrogénio, gás natural e gás de 
petróleo liquefeito. !4. 15) com zeólitos tipo X, 
aproveitando a sua elevada selectividade energé- 
tica de adsorção 

— secagem de etileno e propileno destinados à poli- 
merização, com zeólitos tipo A, aproveitando a 
sua propriedade de peneiro molecular. 


4.1.2 Purificação de gases 


O objectivo é eliminar outras substâncias polares (ani- 
drido carbónico, anidrido sulfuroso, etc.) em muitos ca- 
sos simultaneamente com a secagem. 

Como principais utilizações industriais poderemos re- 
ferir a descarbonação do ar, purificação de gases indus- 
triais contendo vapores sulfurosos e nitrosos e vapor de 
mercúrio, dessulfuração do hidrogénio e L.P.G.0% 15, 
produção de hidrogénio de elevada pureza (99,999%) e 
produção de oxigénio e azoto em que o principal parâ- 
metro operatório é a variação da pressão dos gases. 

Uma aplicação que poderá ser de muito interesse é 
no enriquecimento de oxigénio em aparelhos portáteis 
para doentes respiratórios crónicos!º. 


4.1.3 Separação de hidrocarbonetos 


Podemos ter também aqui dois tipos de separação, 
com base na selectividade geométrica ou, na selectivi- 
dade energética de adsorção. 

Como exemplo do primeiro tipo, pode referir-se a se- 
paração de parafinas lineares das isoparafinas com zeó- 
litos 5 A (processos B.P., Molex U.O.P., Iso Siv Union 
Carbide, etc.). 

Do segundo tipo, há a salientar a separação de ole- 
finas (processo Olex U.O.P.) baseada na maior adsorção 
dos hidrocarbonetos menos saturados e a separação dos 
xilenos (processo Parex U.O.P.), baseada na diferente ba- 
sicidade dos diversos isómeros (o metaxileno é o mais 
básico, o paraxileno o mais ácido). 


4.2 Aplicações como permutadores iónicos 


Apesar de actualmente as grandes aplicações dos zeó- 
litos serem no domínio da química dos processos cata- 
líticos e em adsorção, uma das suas propriedades mais 
importantes que é a capacidade de permuta iónica pode, 
como vimos anteriormente, ter aplicação directa. 

Já na introdução histórica foi referido que a primeira 
aplicação industrial em permuta iónica, consistiu na re- 
cuperação de isótopos radioactivos, tirando partido da 
grande estabilidade dos zeólitos às radiações ionizantes 
e às altas temperaturas resultantes do decaimento de es- 
pécies altamente radioactivas. Neste domínio, deve sali- 
entar-se a grande importância que os zeólitos tiveram na 
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resolução dos problemas de contaminação radioactiva, 
no acidente nuclear de Three Mile Islande. 

No início dos anos 70, outra aplicação surgiu no tra- 
tamento de esgotos para eliminação da amónia, tirando 
partido da elevada afinidade de alguns zeólitos para os 
catiões amónio. 


Uma outra aplicação relativamente recente e em 
grande expansão, é a sua utilização como coadjuvantes 
de detergentes domésticos, para remoção de catiões como 
o cálcio e o magnésio que conferem características de 
dureza à água. Os zeólitos têm vantagens em relação aos 
desendurecedores comerciais à base de compostos orgã- 
nicos, já que não são poluentes e são mais selectivos, 
não tendo maior custo!”, 


Este grande impulso na utilização dos zeólitos nesta 
área, deve-se em grande parte às medidas rigorosas de 
protecção do meio ambientel*, 


A utilização dos zeólitos, como permutadores ióni- 
cos está actualmente numa fase de grande expansão, po- 
dendo referir-se mais algumas aplicações muito promis- 
soras: 


— Imobilização definitiva de compostos radioacti- 
vos de longa vida. Após a sua incorporação no 
zeólito, este seria calcinado a elevada temperatura, 
obtendo-se um material vítreo altamente imper- 
meável. 

A utilização dos zeólitos para retirar amónia dos 
esgotos urbanos, pode ser estendida a esgotos in- 
dustriais e agrícolas. Efectivamente, é possível ob- 
ter uma combinação dos zeólitos com bactérias 
nitrificantes que permitem uma conversão eficiente 
da amónia dos esgotos em nitratos, apresentando 
esta “simbiose” um elevado efeito sinergético. 
A elevada selectividade de alguns zeólitos para a 
amónia permite também a sua utilização em “rins 
artificiais”, possibilitando a fabricação de apare- 
lhos de diálise portáteis": 19, 

Utilização como fertilizantes, dada a possibilidade 
da sua ligação “'simbiótica”” com microrganismos 
nitrificantes e a facilidade de transportarem ca- 
tiões imprescindíveis ao crescimento de plantas 
(macronutrientes como o sódio e o potássio e mi- 
cronutrientes como cobre, ferro ou manganês!?)). 
Utilização na alimentação de animais, melhorando 


as condições de saúde por um efeito aparente- : 


mente antibiótico. Nos animais ruminantes faci- 
lita a absorção de azoto não proteico, técnica já 
anteriormente conhecida, mas de difícil aplicação 
até ao aparecimento dos zeólitos como transpor- 
tadores de iões amónio, 


As aplicações anteriores referiram-se sempre à utili- 
zação dos zeólitos sintéticos. Os zeólitos naturais, par- 
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ticularmente em locais próximos dos jazigos, são utili- | 


zados como suplemento na produção de papel e, 
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também, na indústria da construção civil (Europa), in- 
corporados em cimento!” e em materiais cerâmicos. 


4.3 Aplicações em catálise 


Os zeólitos podem catalisar um elevado número de 
reacções químicas. Como já referimos anteriormente, os 
canais e cavidades dos zeólitos têm dimensões molecu- 
lares, podendo-se assim regular o acesso das moléculas 
aos centros activos. Essa propriedade designa-se por se- 
lectividade de forma, englobando além de efeitos pura- 
mente geométricos, efeitos cinéticos ligados à difusão. 

Os zeólitos possuem elevada área específica e esta- 
bilidade térmica, podendo as suas propriedades quími- 
cas serem modificadas por permuta iónica ou tratamen- 
tos térmicos. Assim, conforme o catião de compensação, 
podem apresentar centros ácidos que catalisam reacções 
de cracking e isomerização, ou centros metálicos catali- 
sadores de reacções de hidrogenação. 

Embora exista, como vimos, uma grande variedade 
de zeólitos naturais e sintéticos só um número muito re- 
duzido teve sucesso como catalisadores, sendo de salien- 
tar o zeólito Y, a mordenite e, mais recentemente, os zeó- 
litos da série ZSM. 

O conjunto das propriedades acima referidas permite 
a aplicação industrial dos zeólitos, como catalisadores 
numa elevada gama de reacções, com saliência para os 
processos catalíticos envolvidos na refinação do petró- 
leo e em petroquímica. 

Descreveremos a seguir as principais aplicações in- 
dustriais em catálise. 


Cracking catalítico 


Constitui a mais importante utilização dos zeólitos 
em catálise, obtendo-se menos coque, menos produtos 
leves e mais destilados médios do que com os catalisa- 
dores amorfos de sílica-alumina, o que conduz a eco- 
nomias consideráveis de energia. 

Os zeólitos utilizados são do tipo Y, adicionando-se 
catiões de terras raras (cério e lantânio) que conferem 
uma maior estabilidade e actividade), 

Aproximadamente 90% das unidades industriais que 
utilizam zeólitos são do tipo FC.C. (fluid catalytic crac- 
king), com tempos de contacto catalisador-hidrocarbo- 
netos de algumas dezenas de segundo que diminui a for- 
mação de coque e evita a degradação dos produtos pelo 
cracking secundário, melhorando-se assim o rendimento 
e a qualidade das gaolinas. 


Hidroisomerização 


A necessidade de obtenção de gasolinas em C, e €, 
de índice de octano elevado, sem adição de compostos 
de alquil-chumbo, exige que a mistura contenha aromá- 
ticos e isoparafinas. A título de exemplo, podemos refe- 
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rir que o índice de octano do n-pentano é 61,7 e o de 
isopentano 92,3. 

Torna-se, pois, necessário em relação às isoparafinas, 
proceder à transformação dos C,/C, que saem da refi- 
naria, isomerizando o n-pentano e o n-hexano!?? 23, 24), 

O processo industrial mais conhecido é o Hysomer 
da Shell) que utiliza catalisadores zeolíticos bifuncio- 
nais, à base de mordenite na forma ácida com 0,5% (em 
peso) de platina. O rendimento pode ser melhorado se 
ao processo catalítico for acoplado uma unidade de ad- 
sorção selectiva com zeólitos que faz a separação das 
i-parafinas das n-parafinas. O conjunto dos dois pro- 
cessos!2%) é denominado processo TIP (total isomeriza- 
tion process) que permite obter uma mistura de hidro- 
carbonetos com um índice de octano de aproximada- 
mente 90. 


Hidrocracking 


Desde a crise petrolífera de 1973, com a necessidade 
de converter uma maior quantidade de produtos petro- 
líferos, várias unidades de hidrocracking têm sido im- 
plantadas na Europa e, em especial, nos Estados Unidos. 

O sucesso deste processo!?? deve-se em grande parte 
à utilização de catalisadores zeolíticos bifuncionais (metal 
nobre-ácido) que apresentam em relação aos catalisadores 
à base de sílica-alumina amorfa, uma maior actividade 
e uma elevada resistência aos compostos azotados. 

Recentemente, para melhorar a resistência à desacti- 
vação pelo coque, têm vindo a ser efectuadas investiga- 
ções” com catalisadores zeolíticos bimetálicos. 


Produção de gasolina a partir do metanol 


Este processo constitui uma das mais recentes apli- 
cações!??) dos zeólitos tipo ZSM, que são muito selecti- 
vos, obtendo-se uma gasolina rica em aromáticos e ole- 
finas, não contendo hidrocarbonetos mais pesados do 
que os tetrametil-benzenos. 

Este tipo de zeólitos devido à sua estrutura particu- 
lar sem grandes cavidades, apresenta uma elevada resis- 
tência à desactivação pelo coque. Efectua, além disso, 
o chamado “'controlo do tráfico molecular”"C9, impe- 
dindo a contradifusão dos reagentes. 


Conclusões 


Os zeólitos são materiais com um futuro extrema- 
mente promissor, pois além das suas actuais aplicações 
em adsorção, permuta iónica e catálise, muitas outras 
aplicações são previsíveis. 

Existe, um número elevado de centros de investiga- 
ção universitários e industriais, trabalhando no sentido 
de melhorar as propriedades dos materiais zeolíticos já 
existentes, por alteração da sua estrutura e composição, 
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e, também, para obtenção de novos materiais capazes de 
satisfazerem às exigências impostas pela evolução tecno- 
lógica dos processos industriais. 
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A TRANSMISSÃO DO CALOR 


EM RELATIVIDADE 
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RESUMO: 


O estudo relativista da transmissão do calor nos sóli- 
dos não pode ser isolado do estudo das suas deformações 
dado em Relatividade todos os sólidos, incluindo os rígi- 
dos, serem deformáveis. (Estão votadas ao fracasso todas 
as tentativas para estabelecer uma equação relativista da 
transmissão do calor em “sólidos indeformáveis”). 

As equações conjuntas e invariantes do movimento e 
da evolução da temperatura numa barra elástica em que 
há transmissão do calor (incógnitas T = T(x,x,); X= X 
(x,x,)) podem ser obtidas a partir das equações: 


3 (Pot DO 4 pº po? 


| a ) + 
ra (PP HH) + Broto 
Ls E 
É 
i a BPo + (1+B?) A + B Pot) 


1-8 


2 
do B+ pos) 


+; (É Poe : 
fá 


que traduzem a conservação da quantidade de movi- 
mento e da energia (X coodenadas dos pontos da barra; 
Po: Po» do: pressão, densidade e fluxo do calor no refe- 
rencial próprio local). 

Subsiste o problema da escolha das relações constitu- 
tivas: 
HI Po = Poís, T); Po=Po(s, 1); K = K(s, T) 

( s deformação) e da escolha entre as equações: 


E 


Original recebido para publicação em 8/6/85. 


ABSTRACT 


Heat transmission in Relativity 


The relativistic study of the transmission of heat in 
solids cannot be isolated from the study of their deforma- 
tions as in Relativity all the solids including the so called 
rigid bodies are deformable. (Any attempt to establish a 
relativistic equation for the transmission of heat in “inde- 
formable solids” is doomed to failure). 

The simultaneous and invariant equations of the move- 
ment and the temperature of a bar with heat transmission 
(T=T(x,x,); X=X(x, x) ) may be obtained from the 
equations 1 and Il which correspond to conservation of 
both the linear momentum and energy. (X denotes the 
coordinates of the material points of the bar; Po, pos do 
denote pressure, density and heat flux in the local coordi- 
nate system). 

The problem of choosing constitutive relations III and 
of choosing between the equations IV remains open. 


A TRANSMISSÃO DO CALOR 
EM RELATIVIDADE 


A transmissão do calor é de todos os grandes capi- 
tulos da Física abordáveis em termos de equações às 
derivadas parciais o único até hoje “rebelde” a um tra- 
tamento relativista. 

A possibilidade de transmitir instantaneamente sinais 
(provocar alterações da temperatura) a distâncias afas- 
tadas, prevista pela equação de Fourier, choca-se com o 
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valor limite c das velocidades de propagação das parti- 
culas e dos sinais imposto pela Relatividade. 

O natural numa situação como esta (a ideia que 
primeiro ocorre) é a de substituir a equação de Fourier 
por outra tida por mais exacta, que satisfaça às exigên- 
cias da Relatividade e com ela se reencontre no limite 
c—0o, Vários autores se empenharam nesta tarefa ao 
que parece sem resultados convincentes. (Informações 
sobre o assunto podem, por exemplo, ser encontradas 
no livro relativamente recente de A. Bressan [1]). 

Note-se que nesta abordagem o problema não é o 
da exigência de uma invariância relativista. 

Sendo (x,t) as coordenadas de um referencial S em 
que esteja imóvel uma barra, (x';t) as coordenadas de 
um outro referencial S em movimento em relação ao 
primeiro com a velocidade v, e T a temperatura medida 
com termómetros fixos na barra, nada exige que as 
funções T(x,t) e T(x',t) verifiquem a mesma equação 
(ou, por outras palavras, que a equação de T(x.t) fique 
invariante na transformação de Lorentz (x.t) — (x';t) ). 

Com efeito, no segundo referencial a barra está em 
movimento com a velocidade -v e é mormal que esta 
velocidade figure na equação de T(x',t). A barra desem- 
penha aqui o papel de “ether” — o tal “ether” que não 
foi encontrado para o Electromagnetismo e cuja não 
existência se traduz pela invariância das equações de 
Maxwell. 

Penso saber a razão do fracasso das tentativas refe- 
ridas e da situação de impasse em que se encontra a 
teoria relativista da transmissão do calor. 

A explicação parece de resto simples: à boa maneira 
de Fourier os relativistas continuaram a estudar a trans- 
missão do calor numa barra indeformável (correspondente 
ao modelo de corpo rígido-indeformável de Born (1909), 
modelo que deu origem a inúmeros paradoxos em Relati- 
vidade). 

Críticas relativamente recentes (ver por exemplo [2] e 
[3]) parece mostrarem, com toda a evidência, que este 
modelo de corpo rígido indeformável deve ser substituído 
por modelos de corpos rígidos deformáveis dotados de 
convenientes leis elásticas. 

No caso da transmissão do calor, sendo a barra 
elástica, ela deve sofrer deformações sempre que haja 
variações do fluxo de calor dado que a elas correspon- 
dem variações na transmissão da energia e da quanti- 
dade de movimento. 

Esta interacção entre transmissão do calor e elasti- 
cidade, que não pode ser ignorada em Relatividade, 
nem no caso limite dos corpos rígidos (e com mais forte 
razão no caso dos corpos menos rígidos), não foi, que o 
saibamos, até hoje considerada. 


W 
W W 


Como escrever então as equações conjuntas da evo- 
lução da temperatura e do movimento de uma barra 
elástica no caso de haver transmissão do calor? 

O problema pode ser resolvido com uma extensão 
do método (ver por exemplo [4]) que permite no caso 
adiabático chegar por duas vias diferentes (a partir da 
conservação da energia e a partir da conservação da 
quantidade de movimento) à equação do movimento de 
uma barra. 

Sejam (x,t) as coordenadas de um referencial S (qual- 
quer) e X as coordenadas “fixas” dos pontos da barra 
deformável. O movimento e a evolução da temperatura 
podem ser descritos por: X= X(x,t)e T = T(x,t). 

Representemos por go (Xt) o vector fluxo de calor 
no referencial próprio S' da secção X da barra no ins- 
tante considerado. O seu significado é o seguinte: cada 
secção X da barra é atravessada no referencial S' que a 
acompanha (e que regra geral muda de instante para 
instante) pela energia: dE = qoA dt, (dtp tempo pró- 
prio de S'; dt; = 1-B? dt; B' = +: Y velocidade de 
S' em relação a S; A secção da barra). O vector qo é um 
fluxo de calor dada a secção X estar imóvel em S' e ser 
portanto nulo (em S') o trabalho transmitido de um 
para outro lado de X. Usando X = x(x,t) podemos escre- 
ver do = Jo(X.t). 

As componentes do tensor impulsão energia em X no 
instante considerado nas coordenadas (x',x4=ct”) asso- 
ciadas ao referencial próprio local S' escrevem-se: 


Te B= Po O O do 
C 
O 6) O O 
O O 6) O 
|) 
do SM gr 
Cc 


Num material sem “memória” po e Po (pressão e 
densidade no referencial próprio S'”) só devem depender 
da temperatura e da deformação que exprimimos por , 
meio de: 


2) 
Ox dx 


OX | EA A fe 
SK => p=liv= 


- e , NE = 
e dt dx 
a alçõõe. 1, 082 Ox 


Escrevemos pois: 
Po = Po(s, É) 
Po = Poís, T) 


3a) 
3b) 
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Para qo a Física Clássica sugere a fórmula: 


4) GOT Rj 


com K eventualmente dependente de se T. 

Em Relatividade o aparecimento da constante c abre 
a possibilidade de escrevermos fórmulas dimensionalmen- 
te correctas do tipo: 


X c 


em que figura a aceleração Y, no referencial próprio e 
em que R é uma grandeza adimensional, eventualmente 
um número puro. Precisamos assim, de conhecer: 

5) K=K(s,T);R=R(s,T). 


Deixemos de momento de lado a discussão destas 
fórmulas 4 e 4' e da procura das relações constituti- 
vas3es. 

As transformações de Lorentz permitem-nos a partir 
das componentes do tensor impulsão energia em Sº conhe- 
cer as componentes em S. Feitos os cálculos obtemos: 


T=<B 


2q q 
pote. B+B'po” (o Brot(I+B)—E +Bpor 


|-Bº 1-Bº 
o oo o 
6) o oo o 
2 Jo 2 2 Jo 2 
Brot(I+B)-T+Êpos 4 Bpot-—c B+ Poe 
|-pº 1-8” 
As duas equações de conservação: 
7) daT=0 ; daTÊ=0 


traduzem-se por: 


2 
7a) ax(Pot TE B+B' pod), (Bpot(1+BYEHBpoc) o 
E: E 


2 
7b) àx(BPo+U+BITE+Bpoe') , (B'po+2p+poe) 
1-pº 1-Bº di 


No caso de termos resolvido o problema da escolha 
entre as equações 4 e 4'e das relações constitutivas 3 e 5, 
introduzindo esses elementos nestas equações obtemos 
um sistema de equações diferenciais de segunda ordem 
em que as incógnitas são as funções x = X (x,x4) € 
T = (x,x,). É o sistema que procuramos. Note-se que as 
equações 7 são invariantes na passagem das coordenadas 
(x,x,) de S para as coordenadas (x*,x*«) de um qualquer 
outro referencial S*. 


W 
Wo o * 


Voltemos um pouco atrás e consideremos um mate- 
rial (que para simplificar suporemos homogénio) não 
condutor ou muito mau condutor (K=0; R=0) o que 
implica ser nulo ou quase nulo o fluxo de calor (q5=0). 

Admitamos uma situação inicial em ques = le 
T = const. em todos os pontos, a partir da qual o mate- 
rial evolui adiabaticamente (q, =0). Em todos os pon- 
tos a temperatura estará directamente relacionada com s 
pela equação de uma adiabatica T = T,u(s). Introduzindo 
esta relação em 3b obtemos: 


8) Po=Pols Tu(S))= Po, ds) 


que nos relaciona po com s sobre a referida adiabatica. 
(Designaremos por pg a densidade no referencial próprio 
no estado inicial). 

A partir de 8 podemos com facilidade relacionar a 
pressão com s sobre a mesma adiabatica. (Basta relacio- 
nar a variação de po com o trabalho recebido numa defor- 
mação). Temos: 


d 
9) Po=-Pogº- 


Em textos anteriores [1] e [5] estudamos os casos 
particulares: 


pre 


p 
2 = + 


10) Po=5" (1 + 1;Po = 


e 
VI-Fº + F's- F's . 


10”) Po = Ps (1-F?)s 


cmo. nO id Fº é $ 
po=piet ia [1 CR! 


Em que F é um coeficiente característico do material 
(O<F<=1) que respeitam a relação 9 entre po e po. típica 
das evoluções adiabaticas, e que correspondem, o primei- 
ro ao caso limite do corpo rígido (em que as ondas de 
choque longitudianais se propagam com a velocidade c) 
e o segundo aos materiais em que as ondas de choque se 
propagam sem dispersão com uma velocidade inferior a c. 


Consideremos agora o caso adiabatico mais geral de 
Po, = Po, (8) ser qualquer, por ora sem outras restrições 
além de ser uma função decrescente (passarmos a escre- 
ver simplesmente Po, , = Po, (S)). 

Sendo qo = O, as equações de conservação escrevem-se: 


o. [e (Po + Poe?) | + à,, [É Po + poe 1. 


1-8” 
e 
Hb a, [Ponte | a, [a focoe | 


ABRIL 87 83 


Usando 9 para eliminar p, obtemos: 


Po dPo do 
I2a R | 94 [Pas | q 
| 2 |-Bº 
dp 
e 7 4Do 
dx| É SPos | 4; Ba sS|=0 
I2b fe ds | “ |-Pot 1-6? 
Desenvolvidos os cálculos obtemos: 
»n Os 2B OB 
13a) a [2-8 dE asd = EM + 
E êr)ts [E tea J-o 
O 1-6 9B ds 
13b) a[8=; LR dx, +28 E ei 


Ee 9 
(para simplificar escrevemos a — Po. ). 
S 
Utilizando agora a expressão 2 de S e usando as 
notações: 


| dx dx Cia 
14) V e dx 5 Va — dx, , Vu “de , 
E sx dx 
Vs=5> ;Va= Vuc 
dx Oxdx, 


obtemos após cálculos triviais mas bastante longos: 


2 
isa é 


15a)- | 


lt 


-V, 


ca Ma 
— 6 Via +] + na 2 Ma A + 
», 
ru] 
LV Vu 
6 WU |+-: afecande + 
V | 
+U yr ]|=0 


Vemos, assim, que a partir da conservação da quanti- 
dade de movimento e a partir da conservação da energia 


chegamos, não exactamente à mesma equação, mas a 
duas equações que diferem nos factores diferentes que 
multiplicam uma parte comum. Note-se que, como não 
podia deixar de ser, todas as configurações estacionárias 
x=x(x); V,x=0O verificam a segunda equação, a da 
conservação da energia, mas regra geral não verificam a 
primeira. (Verificam a primeira as soluções x=x(t); 
V, Xx=0; s=0O, só interpretáveis como barras compri- 
midas ao extremo que se deslocam com velocidade infi- 
nita, não aceitáveis pois na Relatividade corrente). 

Tirando estes dois casos extremos, as soluções de 15a 
e |5b são comuns e verificam a equação que escrita em 
notações correntes toma a forma: 


(SEy (537 
ox Ox, ] dx O xa 
] 2 (ae D+ cd (142) 
16) = 
E dx !- 2 dx 
=. e Ox, , da X E o) X dx 
dx, Ox dx ds dx? dxdx E 
x x 
xy 
dx (PEy 
Er 


com s dado por 2. Para termos a equação correspon- 


Po. 
dente a um caso particular basta precisar a = Er Poa). 


(A equação correspondente ao caso particular 10" foi 
apresentada em [4]). Encontrada uma solução da equa- 
ção 16 a equação de estado 3a permite em seguida deter- 
minar T=T(x,x4). 

Em linguagem algo simplificada podemos dizer que 
encontramos a equação 16 por dois caminhos. De facto, 
o que encontramos foi um sistema de duas equações 
iguais (diferindo embora num factor). Este sistema em 
que só figura X = X(x,x,) transforma-se no caso q,%0O, 
(K£O, RO) no sistema de duas equações diferen- 
tes obtido a partir de 7 em que figuram as incógnitas 
x=x(x;x,))e T= T(x,x,) que só podem, em consequên- 
cia, ser determinadas em conjunto e não isoladamente 
uma da outra. 


*r 
WoW 


No que diz respeito às relações constitutivas 3 e 5 o 
problema põe-se do seguindo modo: 
São admissíveis as relações expressas por quaisquer 
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funções? Há limitações impostas pelo seu significado físico 
(como é, por exemplo, a condição de a função apresen- 
tada em 8 ser decrescente)? Há entre elas relações impos- 
tas por princípios gerais (como a relação 9)? Quais são as 
leis correspondentes aos casos limites e particularmente 
simples (como o são, por exemplo os casos 10 e 10” na 
hipótese adiabática)? 

São problemas em aberto. 

Um outro problema em aberto é o da escolha entre as 
fórmulas 4 e 4”. 

Estando a temperatura relacionada com a energia e 
tendo a energia massa deve sofrer o efeito da aceleração. 
(Em particular, o equilíbrio térmico é o equilíbrio de deter- 
minado tipo de energia que tem massa e como tal deve ser 
influenciado pelas acelerações). 

Uma experiência para averiguar qual das duas fórmu- 
las é a correcta pode ser a seguinte: num disco com ter- 
mómetros em vários pontos basta verificar, atingida a 
situação estacionária (q, = O), se a temperatura é a mesma 
em pontos situados a diferentes distâncias do centro. 
Caso sim é válida a fórmula 4, caso não a fórmula 4" com 
RXO. 

Não cremos que os experimentadores consigam nos 
tempos mais próximos precisões suficientes para detectar 
eventuais diferenças de temperatura... sem partir o disco. 

Restam-nos experiências conceptuais. 

O que é um termómetro? Podemos imaginar que no 
interior de um material com determinado estado de ten- 
são e temperatura é feito um pequeno cavado esférico 
(fig. 1) tendo no seu interior uma barreta ab ligada ao 
material por uma ponte condutora cd. O comprimento 
“da barreta sendo independente do estado de tensão pode 
servir para medir T. Naturalmente que a temperaturas 
mais altas correspondem maiores quantidades de energia 
na barreta. 


NS 


Admitamos agora que numa barra que roda em torno 
de um ponto O há nos pontos A e B(OA < OB) dois 
cavados com barretas. A maiores ou menores temperatu- 


ras correspondem maiores ou menores quantidades de 
energia “térmica” acumulada nas barretas onde não há 
energia elástica. Na situação estacionária o equilíbrio 
térmico mantém-se entre regiões vizinhas. Mas manten- 
do-se tudo na mesma tudo se passa como se o equilíbrio 
se estabelecesse directamente entre a barreta colocada em 
A e a barreta colocada em B. Parece, assim, que nas 
situações estacionárias a distribuição de temperatura pode 
ser analizada independentemente do estado de tensão. 
Ora tendo a energia “térmica” massa as acelerações ten- 
dem a transmiti-la de A para B razão pela qual a tempe- 
ratura em B deve ser superior à temperatura em A (e 
válida a fórmula 4” com R negativo). Na verdade talvez 
seja só a intuição, estado embrionário do raciocínio, o 
que nos faz preveligiar a fórmula 4º relativamente à 4. 

Não pensamos que os matemáticos se possam seria- 
mente ocupar nos tempos mais próximos do estudo de 
equações como a equação 16. No caso a três dimensões a 
situação é muito mais complicada. No caso adiabático, 
postuladas as relações constitutivas e escrito o sistema de 
equações conjunto para T = T (x; x,) e para os três 
XX, (x;, X,), que ainda não sabemos escrever mas que 
alguém escreverá um dia, o estudo das suas soluções 
afigura-se bem difícil. Mas mesmo antes desse trabalho 
estar feito um físico pode desde já ter a convicção de que, 
sob reserva das relações constitutivas terem sido conve- 
nientemente escolhidas, no final, se irá encontrar o limite 
c para a velocidade de propagação dos sinais que os 
relativistas procuram desde o início incluir na teoria da 
transmissão do calor. Não se trata neste caso de uma 
questão de intuição, mas de algo mais sólido, é uma ques- 
tão de confiança, de “fé”, na coerência global das cons- 
truções baseadas nos grandes princípios. 
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VIABILIDADE ECONÔMICA DE BATERIAS 
DE LÍTIO PARA VEÍCULOS ELÉCTRICOS 


C.AC. de Sequeira 


Laboratório de Electroquímica, Instituto Superior Técnico, Lisboa, Portugal 


RESUMO 


A mutação que se perspectiva para a indústria de 
transportes rodoviários com o desenvolvimento de veículos 
eléctricos tem interessado os poderes públicos nacionais 
e comunitários da Europa, os quais têm vindo a investir 
elevadas somas de dinheiro em projectos de investigação 
e desenvolvimento no sector dos veículos de tracção eléc- 
trica. 

No ano lectivo de 1981/82, o autor desta comunica- 
ção esteve integrado num desses projectos. Por essa oca- 
sião iniciaram-se estudos sistemáticos duma célula total- 
mente sólida de Li/polímero/V,O,, bem como dos 
materiais que a constituem. 

O conceito duma bateria de tracção totalmente sólida 
é revolucionário pelo que é prudente fazer uma pausa no 
seu desenvolvimento e aproveitar essa pausa para tentar, 
embora duma forma necessariamente preliminar, anati- 
sar a viabilidade tecnológica e estimar os custos de ma- 
nufactura envolvidos num tal projecto. E essa tentativa 
que se apresenta aqui. Baseados nos resultados labora- 
toriais obtidos com pequenas células iguais à que se es- 
quematiza na Fig. 1, bem como em informação colhida 
junto de fabricantes de produtos químicos e de baterias 
e de industriais do sector automóvel, foi possível conce- 
ber uma bateria que, uma vez devidamente desenvolvida, 
será um concorrente sério face a outras baterias de trac- 
ção e irá brevemente servir de ponta de lança do desen- 
volvimento acelerado do sector de viaturas de tracção eléc- 
ÍrICa. 


ABSTRACT 


In this paper an outline design study has been made 
for an all-solid-state lithium electric vehicle traction bat- 
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tery based upon an extrapolation of the results for small- 
-scale laboratory cells. With a bipolar battery design it 
is shown that very high energy densities (425Wh/kg and 
530 Wh/l) are potentially attainable at realistic-power den- 
sities. Preliminary calculations indicate a materials of 
construction only cost of around 6000 escudos/KWh for 
the battery modules. This figure may be reduced to 4000- 
-5000 escudos/KWh by lowering the cost of the lithium 
salt in the polymer electrolyte and of the intercalation 
compound in the cathode. It is concluded that on pre- 
sent technical and cost considerations the selling price of 
traction batteries for some existing electric vehicles is eco- 
nomically viable, and this novel battery concept holds con- 
siderable promise and merits further research and deve- 
lopment. 


CÉLULA LABORATORIAL MODELO 


O sistema electroquímico ''*mais prometedor”, ou 
seja, que presentemente serve de base ao nosso estudo 
conceptual, é do tipo: | 


V,O; (compósito 
polimérico) 


Electrólito polimérico 
(PEO),LiF,CSO, 


Lítio” 
(ou Li/A1) 


em que o electrólito é uma solução sólida de sulfonato 
de trifluormetano de lítio em óxido de polietileno. O cá- 
todo compósito de óxido de vanádio, V,O,, é composto 
por V,O,, (22,5%), polímero (25,0%) e negro de aceti- 
leno (2,5%); a porosidade quantifica os restantes 50% 
da sua composição volumétrica. A célula descarrega-se 
reversivelmente a 120-140ºC de acordo com a reacção: 


BLi + VOO LisVçO; 
variando as tensões aos terminais entre 3,0-1,7 volts ou 


2,7-1,4 volts, consoante se trate de ânodos de Li ou Li/Al, 
respectivamente. 
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ro positivo 
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Electrólito Anodo | : 
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Fig. 1 Célula laboratorial de Li/polímero/V,O,; 


Com células unitárias consistindo de 


Li (500 um) | Polímero (25 um) V,O (50 um) 

foi possível obter rendimentos catódicos iniciais da or- 
dem de 100% para descargas de 0,25 mA/cmº (duração 
de descarga: 10 h) e da ordem de 50% para descargas 
de 0,8 mA/cm? (duração de descarga: 3 h). Foi também 
possível manter correntes de pico da ordem de 1,5 mA/ 
/em? (regime de descarga: 1,7 h) durante vários mi- 
“nutos. 

Para os objectivos do presente estudo, adoptou-se 
uma densidade de corrente de 0,25 mA/cm? e uma du- 
ração da descarga de 10 h. Trata-se de valores compati- 
veis com ciclos de vida longos em acumuladores bem di- 
mensionados. Para o cálculo das energias de pico 
(energias necessárias quando o veículo eléctrico é sub- 
metido a acelerações rápidas), adoptou-se a d.d.c. de 
1,5 mA/cmº. 


ESPESSURA DOS ELÉCTRODOS 


Recorrendo às leis de Faraday, é fácil mostrar que 
para fornecer uma corrente de 0,25 mA /cm? durante 10 
horas a bateria a conceber basear-se-á em “bolachas” 
delgadas com as espessuras aproximadas: 


Li 24 um 
V,O, 50 um 


Electrólito 25 um 
Célula unitária 99 um 


DESIGN DE BATERIAS E CÉLULAS 


As placas das baterias de acumuladores de chumbo 
poderão atingir áreas da ordem de 20 cm x 20 cm. Selec- 
cionaram-se estas dimensões para os eléctrodos e elec- 
trólito da nossa célula. Por conseguinte, a 0,25 mA/cm”, 
a célula unitária debitará uma corrente de 100mA e apre- 
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sentará uma tensão aos terminais de 2,4V (tensão mé- 
dia para o caso de ânodos de lítio). 

Vinte destes acumuladores ligados em série a debi- 
tar 100mA apresentam uma tensão aos bornes de 48V. 
Trinta destas unidades bipolares, ligadas em paralelo, ori- 
ginam uma “bateria unitária”” de 600 células, que debi- 
ta 3A e apresenta uma tensão aos bornes de 48V. Admi- 
tindo que a descarga se processa continuamente durante 
10 horas, a bateria unitária terá uma capacidade de 30 Ah 
a 48V, ou seja, terá uma energia útil de 1440 Wh. 


Ao tentar conceber os recipientes das células assumiu- 
-se que os bipólos metálicos têm a espessura de | um € 
que os colectores de corrente têm 12 um de espessura 
(filme polimérico com 10 um de espessura sustentando 
uma camada metalizada com 2 um de espessura). A es- 
pessura da bateria unitária pode agora calcular-se: 


600 células com espessura unitária de 99 um ...... 5,94 em 
540 bipólos com espessura unitária de 1 um ....... 0,054 em 
60 colectores de corrente c/ espessura 
unitária de E um cursa onrserenenirenapsaseenánia 0,072 em 
Espessura total da bateria........cseesceeecas 6,066 cm 


O recipiente da bateria unitária, fabricado de folha 
de aço inox com 0,25 mm de espessura, terá as dimen- 
sões internas 20,1 cm de largura x 22,5 cm de altura x 
x 6,1 cm de espessura. Propõe-se deixar 2,5 cm de altura 
acima das placas para acomodar as barras de união e 
os bornes. 


ENERGIAS ESPECÍFICAS 


A bateria unitária tem um volume de 2,7 litros 
(20,1 cm x 22,5 cm x 6,1 cm) e uma energia de 1440 Wh; 
por conseguinte, a sua energia específica em volume é 
de 533 Wh/litro. 

Para calcular a energia específica em peso há que con- 
siderar a massa dos constituintes individuais da bateria 
unitária. Estas são como se indica a seguir: 


E a ARES MIRA 1 (7789, ON 304 g 
Maid poi e EA PRP OSC AV 720 e 
Compósito de VOO sessssssisaasasisios ails Goa de inata ças 1946 g 
POLOS MEANCOS usas penas 53 


SECR GM DAN qa se SS SS 293 p 
Barras de união de CODre.....eesmenesacespeseseesereeasesso 32 8 
FEMALE: DO CODIE usasenceeerssatesarterios presas arts pepms 25 e 

Prso-total da BALOTia essas asas sasasaarace 3378 & 


A energia específica em peso é, então, igual a 1440/ 
/3,378 = 426 Wh/kg. 
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POTÊNCIAS ESPECÍFICAS 


A potência debitada continuamente (duração de des- 
carga: 10h) pela bateria é de 144W (3A e 48V). Con- 
tudo, a potência de pico debitada (assumindo-se a d.d.c. 
de 1,5mA/cm?) é de 864W para acelerações rápidas 
dos veículos eléctricos. 

Um parâmetro a ter em consideração no design de 
uma bateria de tracção eléctrica é a razão energia/po- 
tência de pico. Esta razão deve ter um valor de cerca de 
| para uma furgoneta urbana e de cerca de 3 para um 
autocarro urbano. Para a nossa bateria, esta razão é 1,66 
o que se ajusta perfeitamente com o objectivo para que 
se concebe. 


BATERIA PARA VEÍCULOS DE TRACÇÃO 
ELÉCTRICA 


Existem quatro viaturas europeias que são comercia- 
lizadas em duas versões: a de veículos eléctricos e a de 
veículos de motor a gasolina ou Diesel. Presentemente, 
as versões eléctricas são economicamente desfavoráveis 
face às de combustão interna, mas esta situação alterar- 
-se-á assim que for possível aumentar o tempo de vida 
das baterias ou reduzir o seu preço de venda. A substi- 
tuição da bateria de Pb/ácido pelo novo conceito de ba- 
teria aqui delineado constituirá um primeiro passo nesse 
sentido. Para três dos veículos referidos, o uso da nova 
bateria levará às seguintes situações: 


Veiculo eléctrico Vw GOLF PH-FIAT HH E DAIMLER-BENZ 207 
Peso das baterias de 320 ke 400) kg I000 kg 
Pb/ácido usadas 

lensão nominal 96 V 84 V I8O W 


Baterias unitárias 2 pe 4 
(192V) 


ligadas em série (96V) (96V) 
Capacidade nominal 130 Ah 160 Ah 160 Ah 
Baterias unitárias 4-5 á Í 
ligadas em paralelo (120-1504h) (150 Ah) (150Ah) 
Nº total de baterias S-10 LO 20 
Peso das baterias de 27-34 kg 4 ke 68 kg 


Li propostas 


Trata-se de configurações sugeridas, a titulo ilustra- 
tivo, mas está claro, que existem outras possibilidades. 
Por exemplo, a tensão aos terminais da bateria pode 
ajustar-se por modificações do número de acumulado- 
res ligados em série e em paralelo. 


CUSTOS DOS MATERIAIS 


Os custos dos materiais foram avaliados duma ma- 
niera preliminar para o design de bateria descrito acima. 
As cotações globais foram obtidas de fabricantes de pro- 
dutos químicos e as mais favoráveis indicam-se a seguir. 

Estes valores referem-se a custos minimos que não 
entram em conta com desperdícios e sobras (p/sucata). 
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A eles há que acrescentar os custos de fabrico e de mon- 
tagem; parece razoável considerar que esses custos não 
ultrapassariam os primeiros, de modo que o preço final 
da bateria rondaria os 12 000 escudos/KWH. Após adi- 
ção dos gastos gerais e lucros, o preço de venda ao pú- 
blico deve estar ainda dentro 


Custos das Matérias-Primas 


Matérias-primas g/KWh Esc/KWh 
Lítio metálico 2 OO 
Polimero (PEO) 540 620 
Sulfonato de Trifluormetano de lítio 218 1740 
TO 975 2320 
Negro de fumo 96 60 
Bipólos metálicos 37 20) 
Recipiente da bateria (aço inox) 203 80 
Barras de união e terminais de cobre 40 u 
5960 


da gama de preços praticada com a bateria de acumu- 
ladores de chumbo (12 000-30 000 escudos/KWh). Além 
disso, é admissível a substituição de certos materiais e 
a redução dos seus preços; por exemplo, devem conseguir- 
-se compostos de inserção de lítio mais baratos e tam- 
bém é altamente viável o uso de sais de lítio menos dis- 
pendiosos que o sulfonato de trifluormetano de lítio. 
Quanto às tecnologias de fabrico, deve também ser pos- 
sível usar processos economicamente eficazes, nomea- 
damente no que respeita à obtenção de películas delga- 
das de lítio, de modo a reduzir o preço final da bateria. 

Em conclusão, esta bateria sólida de lítio, que opera 
à temperatura de 120º€, com uma energia especifica em 
volume de 530Wh/1, uma energia específica em peso de 
425Wh/kg, uma razão energia/potência de pico de 1,6- 


-1,7 e preço competitivo deve ser um concorrente sério - 


para o mercado de veículos rodoviários de tracção eléc- 
trica. 
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ON THE CALCULATION OF THE FOUR-VECTOR 
(k, w/c) IN A LINEAR TRANSFORMATION 
OF COORDINATES. AN EXAMPLE. 


Maria Hermínia Marçal* 


ABSTRACT 


The postulate of the general covariance of Maxwell 
equations provides a new insight on the problem of coor- 
dinate transformation. Thus, problems associated with 
inertial frames, such as the propagation, reflection or 
transmission of electromagnetic waves, do not necessari- 
ly require the use of the Lorentz transformation, and it 
can even be demonstrated that for linear transformations 
the formalism of the plane monochromatic wave is still 
possible. 

Although that is not our aim in this article, the study 
of accelerated systems is also based on the general covar- 
iance of Maxwell equations. However, if non-linear trans- 
formations are used, the formalism of the plane mono- 
chromatic wave ceases to be adequate. 


RESUMO 


Postulando a covariância das equações de Maxwell 
abrem-se novos horizontes sobre o problema da transfor- 
mação de coordenadas. Assim, os problemas ligados a 
referenciais de inércia, como seja a propagação, reflexão 
e transmissão de ondas, não exigem necessariamente a 
conservação da transformação de Lorentz, e até se de- 
monstra que para transformações lineares o formalismo 
de onda plana e monocromática ainda é possível. Embora 
isso não seja feito aqui, a covariância das equações de 
Maxwell permite também tratar sistemas acelerados. Mas 
intervindo relações não lineares na transformação de coor- 
denadas, o formalismo de onda plana e monocromática 
já não é possível. 


|. INTRODUCTION 


This article is based on material contained in the dis- 
sertation submitted in partial fulfillment of the require- 
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ments for a Master of Science degree!!!. However, some 
formal passages have been omitted and the sequence and 
critical analysis of the results is different, clearly demon- 
strating how powerful is the postulate of the general 
covariance of Maxwell equations. 

It is indeed an immediate consequence of postulat- 
ing the general covariance of Maxwell equations that a 
plane monochromatic wave maintains its characteristics 
under a linear transformation of coordinates. In fact, 
given two coordinate systems S and S, it results from the 
tensor transformation law together with the uniqueness 
of a linear system's solution that!? 


1) To the amplitude (F,.), of a plane monochro- 
matic wave in S we can associate a certain quan- 
tity (F,), through the relation” 


ox” 
Ox* 


Oxm 


(Edo o axi 


É À (1.1) 


(F,), which is independent of the space-time 
point (x') for linear transformations of coordina- 
tes, will be considered from now on as the ampli- 
tude of the wave in S. 

11) The phase is an invariant, q = é, which, consider- 
ing the tensorial nature of the four-vector (x'), 
defines a new four-vector (k') as follows 

d=kix=kix=ó (1.2) 
with k, and k, related by the tensor transforma- 
tion law. 


However it should be noted that if S is not a carte- 
sian coordinate system, a careful interpretation of the 


(*) Assistente do IST, Dep. Eng. Electrotécnica, Centro de Electrodinâmica 
da UTL, INIC. 

(**) The conclusion is identical for contra-or covariant components. 

Original recebido para publicação em 21/11/85. 
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four-vector (k') is needed. Naturally, the same remark 
can be made about space-time, or any other physical 
quantity. 

Based on the fact that the characteristics of a plane 
monochromatic wave are maintained under a linear trans- 
formation of coordinates, we shall present a systematic 
method for solving the problem of reflection and trans- 
mission of electromagnetic waves by a moving medium. 

As it is well known, when the medium is moving with 
constant velocity the problem is solved using the Lorentz 
transformation between two inertial frames: the labora- 
tory frame | and the moving frame I. When the trans- 
formations involved are linear, though not necessarily 
Lorentz transformations, the moving frame is still iner- 
tial but the coordinate system S associated with | is no 
longer a cartesian system. Therefore, the study of elec- 
tromagnetic phenomena in S will require the postulate of 
the general covariance of Maxwell equations. Neverthe- 
less, because the characteristics of a plane monochro- 
matic wave are preserved, the classic formalism may still 
be used. 


2. THE GALILEAN TRANSFORMATION AND THE 
REFLECTION AND TRANSMISSION 
OF ELECTROMAGNETIC WAVES BY 
A UNIFORMLY MOVING DIELECTRIC MEDIUM 


Let S(X, Y, Z, €T) be the coordinate system associat- 
ed with an inertial frame Iº, We shall call n,, the met- 
ric tensor in S defined by the fundamental quadratic form 


ds? = -dX?-dY?-dZ” + cidTº = mn, dX' dX: (2.1) 
The Galilean transformation 


Y=y 
1.:=1 


X = 1: Vit 
di? 


(2.2) 


is a linear transformation relating the coordinate systems 
S and S. Each reference point of S is seen in S moving 
with uniform velocity v in the X direction. Thus we can 
associate S with an inertial frame 1. 

The metric tensor in S, g,, obtained from relations 
(2.2) either by applying the tensor transformation law 
to n, Or by imposing the invariance of the quadratic 
form ds”, takes the following form 

81=Bp=Bn=-| gu=1-B* B=v/c 

(2.3) 
Z;=2n=-B g,=0, 1%)=1, 2, 3 

The Galilean transformation will be used to solve the 
problem of the reflection and transmission of electromag- 
netic waves by a lossless, non-dispersive dielectric medi- 
um moving with constant velocity. 
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In a previous work!!, the method was presented in 
detail, including calculations for a horizontally-polarized 
incident wave. In order to avoid fastidious repetition only 
the main steps and a summary of results will be presented 
here. 

The geometry of the problem is shown in Fig. 1, and 
O is the incidence angle. In S the incident monochro- 
matic plane wave may be written as follows 


Y X 


Fig. 1 Geometry of the problem 
E, = Ely eis! 


6!=K! Xi=m, K! Xi= (2.4) 


| 
=" aT-(E cosm)x-(- sin 6)Y 
c c c 
In the moving system S, defined by (2.2), the inci- 
dent wave is also a monochromatic plane wave which 
can be expressed as 
E, = El, dé 
(2.5) 
dl=g, ki xi 


with g, given by (2.3) and K! and k! subjected to the 
tensor transformation law. 

The reflected and transmitted waves are then calcu- 
lated following the classical method. Thus, by imposing 
the boundary conditions at the interface between the two 
media, which is at rest in S, we obtain 


Wii) The continuity of the phases of the incident, 
reflected and transmitted waves at the interface, 


B; kli X=g, kRi XM=g, kTi x) 
x! = 0 


iv) The equality of the amplitudes of the tangential 
electric and magnetic fields at the interface, 


(E, + Ed = (Ed 


(H4 + HB, = (HD, 


(2.7) 


(*) Capital letters indicating a (pseudo-) cartesian system. 
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| 


(2.6) | 
| 


Table 1 Inc. Wave 
Disp. Eq. E; ki ki=0 


(a). ER 


Refl. Wave | 


8 kk Kti= =U 


Trans. Wave 


É Eis (º) 


(2º) (b) . EI; (90) 


(92 (1-82 02) TP A KIA (KT +28 nº KT KT (1-02) nº (KT =0 


with nº=eu/€, Ho 2 = (1-52)! 
(**) Q= [n? (1-8 cos 02-(1-62) sin? 0]12 


(*) R,=Iícos 9-8)-QJ/I(cos 6-8) +Q]; T,=1+R, 


(20) (a)=(kRi+ 8 kRº)/[(1-82) kR4-8 KRI]; 


MORTA 
eso 


- To calculate all the components of the four-vectors 
kRi and kf, we use relation iii) together with the disper- 
sion equation in S for each medium. 

The reflection and transmission coefficients, and con- 
sequently the amplitudes of the electric and magnetic 
fields of the reflected and transmitted waves, can be cal- 
culated from relation iv), Maxwell equations and the con- 
stitutive relations of each medium. The results thus 
obtained are presented in Table 1. 

The inverse transformations of (2.2) applied to the 
fields in S finally lead to the reflected and transmitted 
waves in S, as given in Table II. 


3. ANALYSIS OF THE RESULTS 


The results of Table II are well known, today!'!: they 
reflect the changes in the amplitude, frequency — the 
Dóppler effect — and wave vector that in the laborato- 
ry frame are to be expected in the reflected and trans- 
mitted waves by a lossless, non-dispersive dielectric medi- 
um in uniform motion. The same results would have been 
obtained if, keeping within the scope of special relativi- 
ty, we had used the Lorentz transformation*!, instead 
of the Galilean transformation. 
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y2º o! (Il - 2Bcos 6 + 82)/c 


(b) = (KT + 8 kT9/((1-82) kT-8 KT] 


-wl sin 0/c 


ywl (1-5 cos 0) + 8 Q]/c 


Table I presents the results obtained in the moving 
system S (x, y, z, ct). Naturally in another coordinate 
system S' (X, Y, Z, cT), even if moving with the 
same velocity v, the conclusions would be different, and 
the comparison is interesting. 

Let us therefore consider this new coordinate system 
S* and its relations with S 


(3.1) 


It is obvious that the time intervals measured in each 
system at a given reference point (dX =0=dx = 0) are 
related by 


d (3.2) 


Consequently there is no agreement in what concerns 
physical entities depending on time even if — and this 


(*) Again a pseudo -cartesian coordinate system. 
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is the simplest case — the measurements are made at 
the same point. 

It is interesting to note that identical conclusions can- 
not be so easily obtained for spatial measurements, since 
simultaneous events in S' are generally not simultane- 
ous in S. 

In fact, the absence of time-orthogonality of the met- 
ricin S forces us to a re-definition of the spatial distance, 
which will he given by! [5] 


-dl2=(g, -B 8/2) dx' dx) (3.3) 


lj=1,2, 3 
As an example the following situations may be con- 


sidered: 


v) In'S the frequency of the incident wave is (from 
(2.5) and Table 1) 


E =-w! [(cos 0-8) B-(1-B)] = 
=wl (1-8 cos 0) (3.4) 
On the other hand in S' we have!!! 
dd ua. DO! 
ai do (1-8 cos 0) = Fe (3.5) 


The relation above was expected, considering 
(3.2) 

vi) The velocity of an electromagnetic wave 
propagating in vacuum along the positive Y axis, 
measured in S” is obviously c = dY/dT, while in 


S we have! 
dl dy dY 2 
— — e O e —- 3.6 
dd dt df * V 1-8 (6) 


vn) The velocity of an electromagnetic wave propa- 
gating in vacuum along the positive X axis is 
c=dX/dT in S”, while in S it will be 


Amb 
NÉ + B 
with dx/dt obtained from (3.1) and equal to 
c(1-6). 


dio [»2 (dx)2)172 E 


3.7 
dt dt a) 


Thus, an observer in S (x, y, Z, ct), even in the sim- 
ple case of vacuum, finds discrepancies between what 
he measures and what he expects from the Principle of 
Special Relativity. 
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4. CONCLUSIONS 


The method here presented, based on the postulate 
of the general covariance of Maxwell equations, is ade- 
quate for solving problems involving non-cartesian coor- 
dinates. The fact that linear transformations are used 
proves to be an advantage, since it is possible in this case 
to follow the well-known formalism of the plane mono- 
chromatic wave. 

It should be noted that contrary to what happens in 
a cartesian coordinate system, the results thus obtained 
cannot be given an immediate physical meaning. This, 
we think, is clearly evident from the examples given. 

On the other hand, although the plane wave formal- 
ism has to be dropped, the method is easily extended to 
non-linear transformations of coordinates. Consequently 
the study of accelerated systems is possible, within the 
postulate of the general covariance of Maxwell equations. 
This has been done int! for a very simple case of 
“constant” linear acceleration — hyperbolic motion. 
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FORMATION OF DIFFUSION COATINGS ON IRON 
AND STEEL. 
PART I: FORMATION PROCESSES 


C. A.C. Sequeira and C.M.G.8S. Nunes* 
Laboratório de Electroquímica, Instituto Superior Técnico, Lisboa 


ABSTRACT shape and size of the article to be protected, the man- 
ner of fabrication, the subsequent service conditions, etc. 
This paper is primarily concerned with metallic diffu- Diffusion coatings are i) alloy coatings produced at 


sion coatings for the protection of iron and steel; special high temperatures by the inward diffusion of the coat- 
reference is made to the physico-chemical factors deter- ing material into the base metal, and, or 1) composite 
mining which metals are able to form diffusion coatings, electrodeposited coatings which are subsequently inter- 


the methods of application available, and the nature of diffused by thermal treatments. However, the terminol- 
the alloy coatings thus formed. ogy “'diffusion coating”” has never been strictly de- 


finedD. The word diffusion may be defined as 1) the 
spreading of a constituent in a gas, liquid or solid, tend- 
RESUMO ing to make the composition of all parts uniform; or 11) 
the spontaneous movement of atoms or molecules to new 


Estudam-se os recobrimentos por difusão usados na sites within a material; or, following the Shorter Oxford 
protecção do ferro e do aço. Nesta 1.º parte, ocupamo- English Dictionary, Hi) the phenomenon occurring when 
nos dos requisitos físico-químicos necessários para a sua atoms of the same, or different, materials “intermingle””. 
formação, das técnicas de aplicação mais usadas e da This effect, when considering the majority of metals, is 
natureza dos recobrimentos resultantes. only significant at relatively high temperatures (>300º€). 

Thus, diffusion coatings must pertain to a high-tempe- 
rature process. The authors therefore consider that all 
1. INTRODUCTION metallic coatings formed at, or subsequently intention- 
ally subjected to, a high temperature whereby an alloy- 

The application of selected protective coatings ena- layer is formed, constitute a diffusion coating. This defi- 
bles iron and steel components to be used in many, other- nition includes, for example, hot-dipping and cladding 
wise prohibitively aggressive, environments. In this way as methods of forming diffusion coatings, although in 


“the obvious advantages of these cheap unalloyed ferrous each case a significant proportion of the coating layer 


metals (for example, good mechanical properties, ease may be comprised of the relatively pure coating metal. 
of fabrication, etc.) may be economically combined with However the use of a high temperature for the applica- 
the increased resistance, to corrosion heat and/or wear, tion of a metallic coating cannot be used as a defini- 
endowed by such surface treatments. tion in itself; the two-stage process of spray-fusing, par- 
Protective coatings may be classified as either metallic ticularly when applying refractory metals, cannot always 
or non-metallic; the latter class can be further subdivided be said to result in a diffusion coating. 
(for instance, organic and inorganic). The method of The commonly available methods for the application 
coating application is, in many cases, dependent on of metallic coatings can be divided into two categories 


several factors; not only the type of coating but also the (a) low-temperature processes and (b) high-temperature 
processes, although there are certain two-stage techniques 


* Present address: Secção Autónoma de Ciência dos Materiais, Univer- that combine both. Electrodeposition and metal Spráy- 
sidade Nova de Lisboa, Monte da Caparica. ing are examples of low-temperature methods. Diffusion 
Original recebido para publicação em 26/12/85. coatings (for example, as produced by hot-dippine, 


ABRIL 87 93 


vapour deposition, cementation, etc.) afford examples of 
coatings produced by the latter. Their outstanding charac- 
teristic is the metallurgical bond produced between the 
coating and the substrate by diffusion-alloying. 

This paper discusses the physico-chemical require- 
ments for the formation of diffusion coatings and the 
techniques used for that formation. Moreover, a brief 
analysis of the nature of the coatings, irrespective of the 
coating technique used, is presented. 


2. REQUIREMENTS FOR THE FORMATION 
OF DIFFUSION COATINGS 


Gorbunov? has summarized the conditions that are 
necessary to allow the formation of metallic diffusion 
coatings on iron. They are: 


|) The atomic diameter of the coating metal must 
not exceed that of iron by more than 15-16 per 
cent. 
ii) the coating metal must be soluble in iron at room 
and at higher temperatures. 
Wi) A metallic contact must be obtained between the 
substrate and the coating metal. 


“Metallic contact”” is not defined and as diffusion 
can take place between all phases the term elemental con- 


TABLE 1 


Elements (metallic) that are known to form diffusion 
coatings on iron 


Solubility (wtYo) 
Size-Factor 


Aluminium 
Antimony 
Arsenic 
Bervllium 
Chromium 
Cobalt 
Copper 
Gold 
Manganese 
Molybdenum 
Nickel 
Niobium 
Rhenium 
Tantalum 
Tin 
Titanium 
Tungsten 
Vanadium 
Zinc 
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TABLE 2 


Elements (metallic) that do not form diffusion coatings 
on iron (they are all insoluble in iron) 


Bismuth 
Cadmium 
Caesium 
Calcium 
Lead 
Lithium 


Magnesium 
Mercury 
Potassium 
Rubidium 
Silver 
Sodium 
Strontium 
Thallium 


TABLE 3 


Elements (metallic) where formation of diffusion 
coatings on iron is uncertain or unknown 


Cerium 
Gallium 
Germanium 
Hafnium 
Indium 
Iridium 
Osmium 
Palladium 
Platinum 
Ruthenium 
Scandium 
Uranium, 
Yttrium 
Zirconium 


tact would seem to be more appropriate. A favourable ! 


size-factor ensures that diffusion is physically possible, 
but this in itself is notsufficient, even in the presence of 
a concentration gradient at the coating-substrate inter- 
face. 

The driving force for diffusion is provided by the free 
energy of solution of the system. There are anomalies 
to the ““size-factor rule”” (c.f. size-factor vs. solubility) 
and as Rhines*) states it is probable that a diffusion 
process can only affect any important changes in the 
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physical characteristics of the bodies involved (ie. a diffu- 
sion coating can be formed) where there exist, in the cor- 
responding alloy system, either relatively broad ranges 
of solid solubility, or compounds, or both. A summary 
of information that relates to the formation of diffusion 
coatings on iron and steel is given in Tables 1-3. 

Generally, the solubility of elements in iron increases 
with temperature, and this factor, together with the in- 
creasing equilibrium-defect concentration, ensures that 
the diffusivity also increases with temperature. The ef- 
fect of temperature on solubility can also be compared 
with the difference of the coefficient of thermal expan- 
sion of the elements; a large value (e.g. aluminium) in- 
creases the size-factor with increase of temperature and 
can maintain the solubility essentially constant. 

As the formation of diffusion coatings is usually car- 
ried out at high temperatures the relative values of the 
expansion coefficient can also determine the residual 
stress inherent in the surface of the coated article”. 


3. METHODS OF FORMATION OF DIFFUSION 
COATINGS 


The formation of a diffusion coating requires elemen- 
tal contact; such contact points are inhibited by oxide 
films, inclusions in the surface, surface debris, etc. 
Thorough cleaning of the surfaces to be coated is there- 
fore essential and fluxes (or reducing agents) are often 
used for the removal of any surface contamination that 
may form before, or during, coating. 

The surface profile can also affect the formation of 
diffusion coatings in several ways dependent on the 
method of application used. For example, a rough sur- 
face will decrease solid-solid diffusion; however, in solid- 
liquid diffusion elemental contact is assisted and gene- 
rally a thicker, albeit less uniform, coating is formed. Li- 
kewise a surface which has been strained or worked (i.e. 
contains residual surface energy) reacts more vigorously. 
The surface preparation before formation of the coating 
can thus be seen to be an important factor. 

It is convenient to discuss the techniques used for the 
formation of metallic diffusion coatings under three clas- 
sifications: 


1) Solid-solid diffusion. 
1) Diffusion between solid and liquid. 
11) Diffusion between solid and gas. 


This classification does not define the mode of diffu- 
sion that occurs during the formation of the coating; for 
example, in hot-dip galvanizing the growth rate after the 
initial formation of the alloy-layer is dependant on diffu- 
sion in the solid state(”, 


3.1 Solid-Solid Processes 


This class includes cementation processes such as 
chromizing and sherardizing, cladding, and the anneal- 
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ing of electroplates, hot-dipped coatings, slurry or 
sprayed-metal coatings!”, There is a great deal of work 
showing that certain metals can interdiffuse at tempera- 
tures below their melting points merely upon contact” 
2) (c.f. friction and pressure welding) and the impor- 
tance of intimate contact between the two metals has 
been demonstrated. Rigg!) found that diffusion of 
Zn & Fe, at temperatures below the melting point, oc- 
curred only at points where contact had been maintained. 
Masing!!?) showed that increased pressure assisted alloy- 
ing by bringing the metals into more intimate contact. 
Metallic compound formation obviously plays a role, and 
a certain amount of diffusion via non-equilibrium defects 
would also be expected in some cases particularly where 
there is a wide separation between the melting-point and 
the temperature of diffusion. In most cases, the cement- 
ing metal also has an appreciable vapour pressure at the 
process temperatures used, and it is probable that some 
metal transfer also occurs through the gaseous phase. 

Cementation consists of heating the metal to be 
coated while in contact with the coating metal, usually 
in powdered form, to a temperature below the melting 
point of the more fusible of the two. An inert filler 
material (e.g. alumina) is added to the cementation mix 
to prevent any fusing of the powder; this allows higher 
operational temperatures and hence lower process times. 
The cementation process can be carried out in relatively 
simple plant but there is obviously a size limitation, and 
economically a large surface-area to weight ratio is desira- 
ble. Coatings of the higher melting-point metals can be 
obtained by cementation processes using “carrier 
formers” (activators), usually halide salts; these are more 
often used nowadays as they allow lower process tem- 
peratures. These techniques combine characteristics of 
both solid-solid and solid-gas methods and are discussed 
under the latter classification in this paper. 

The other methods in this group are basically simi- 
lar to the cementation process but involve, at least, two 
steps: 


1) Contact of the coating metal with the base metal. 

11) Diffusion, by heat treatment (ie. annealing) to ob- 
tain adhesion, and in some cases cohesion, by al- 
loving”, 


Electrodepositiont!S!), metal-spraying?27, and 
physical vapour depositiont??) (see also Solid-Gas 
Processes) are coating methods in their own right. 
However heat treatment of such coatings is often used 
to obtain better adhesion?? or different surface proper- 
ties09, Slurry, or powder, coatings!39 [c.f. application 
of slip coat during enamellingº”] and electrophoretic 
depositionts-“9) [sometimes supplemented by hydropuls- 
ing densification!*!. *)] are also used to some extent in 
order to apply a controlled, uniform, thickness of coat- 
ing metal to the substrate; these have been developed bas- 
ically for application to metal strip. The whole is then 
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compacted and heat treated under controlled conditions 
to obtain both adhesion and cohesion. In this respect 
Bullough'* has presented singular results pertaining to 
the bonding of aluminium foil and powder to steel strip 
in that the bond was shown to be due to the formation 
of oxides rather than to intermetallic compounds (c.f. 
deposition of Cr/Cr,O, coatings for tin-free steel). 
Metal cladding'* is similar to the process of veneer- 
ing. Metal sheets can be successfully bonded together 
under the simultaneous influence of temperature and 
pressure. Alternatively a duplex ingot can be produced 
which is subsequently rolled or drawn down into the re- 
quired form. The process is limited to those metals of 
similar rolling characteristics and the thickness of the 
coating is governed by the relative amounts of each me- 
tal and the amount of reduction received during process- 
ing. Advantages of this process are the excellent density 
obtained and the possibility of applying thick coatings. 


3.2 Solid-Liquid Processes 


The hot-dipping process is one of the oldest methods 
of applying metallic coatings to other metals. Essentially, 
it involves Immersion of the metal to be coated in a bath 
of the molten coating metal. It is necessary that the two 
metals alloy with each other for the successful produc- 
tion of a hot-dipped coating. In this manner, complete 
initial wetting of the surface of the base metal is ensured 
which maintains uniformity of the coating; furthermore, 
as a result of alloying, the adhering film of molten me- 
tal solidifies without showing any tendency to coalesce 
into globules which would give rise to undersirable dis- 
continuities in the coating. In most commercial processes 
(e.g. coating of aluminium, tin, zinc) alloying proceeds 
readily, but although there is no alloying between lead 
and iron, hot-dipped coatings can be produced by ad- 
ding a mutual alloying element (e.g. tin) either to the bath 
or to the surface of the substrate metal. 

The formation of an alloy-layer on the immersed sur- 
face begins as soon as the surface reaches a temperature 
high enough to allow diffusion to occur; this tempera- 
ture may be below the liquidus temperature. The rate of 
formation of the alloy-layer under a given set of condi- 
tions (e.g. alloying elements, surface profile, heat capac- 
ity of the substrate) is a function of the time and tem- 
perature of immersion. 

On removal from the hot-dipping bath a film of li- 
quid metal is withdrawn on the alloyed surface. 
Hrbek“º has related the physical parameters of the 
process, with the coating thickness, by 


E Bali 


3 Y De 
where h is the thickness (um) 
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W is the withdrawal rate (cm. sec!) 
D is the density 
n is th viscosity 


This outer layer, of essentially pure metal, can be 
reduced in thickness, while still liquid, by various 
methods depending on the type of article. 

Hot-dipping is technologically (and economically) 
limited to metals of relatively low melting point, for in- 
stance zinc, tin, and lead alloys, and to a certain extent 
aluminium coatings. The equipment required for the 
process is mainly dependent on the form of article to 
be coated (i.e. sheet, wire, etc.); it is however also affected 
by the type of coating required. Solid-liquid processes 
also involve dipless metallizing coatingst”. 


3.3 Solid-Gas Processes 


Coatings produced by the reaction of the solid me- 
tal substrate with a gas can be considered in two 


1) Physical vapour deposition [e.g. thermal evapora- 
tiont*. *9, cathode sputteringf”]. 
1) Chemical vapour deposition6”. 


The temperatures commonly employed in physical 
vapour deposition generally preclude the possibility of 
extensive solid-state diffusiont*: *9, 

The rections of chemical vapour deposition can be 
effected in several ways: 


|) Thermal decomposition or pyrolysistt, 
1) Hydrogen reduction. 
11) Reduction with metal vapour. 
iv) Reaction with the substrate. 


The substrate may also enter into the reaction in var- 
IoUS ways: 


1) Acting as a reducing agent for a metal halide. 
1) Acting catalytically on the deposition reaction. 
1) Dissolving the coating metal to form a solid so- 

lution or compound. 


Coatings formed by chemical vapour deposition are 
produced by the pack, fluidized-bed!%!. 2), and vapour- 
streaming methods. Pack cementation techniques have 
been discussed extensively by Aves and Ecord'*». The 
basic metal is packed into a powdered mixture composed 
of the metal to be diffused, an inert filler, and an acti- 
vator (usually a halide salt). The use of an activator per- 
mits a lower temperature to be used than would other- 
wise be required by a conventional cementation process. 
The sealed pack is then held at a given temperature in 
a furnace for a definite time. Parts to be coated by the 
fluidized-bed method are suspended in a bed of the coat- 
ing metal which is fluidized by a halogen-gas stream 
mixed with an inert carrier gas. This is similar to the 
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flowing gaseous methods where the gaseous metal halide 
(of the coating metal) is streamed over the metal to be 
coated. 

Salt bathst** have also been used for the formation 
of diffusion coatings; these contain the coating-metal 
halide together with a salt mix which is fused at the 
process temperature. Although allowing extreme opera- 
tional flexibility a serious drawback to its use 1s the ther- 
mal shock received by the component to be coated/?), 
This could cause distortion on entry to the bath, and 
spalling of the coating on with-drawal. This method has 
been extended by the ““electrolytic metalliding” 
process'%) utilising molten fluorides; deposition, and 
diffusion-alloying, occur at the cathode. 

Most of the methods of application of metallic coat- 
ings are summarized in Table 4. The method used has 
a marked influence on the nature of the coating; its thick- 
ness, composition, uniformity, continuity, etc., are all af- 
fected, likewise protective quality, appearance and cost. 

These aspects, particularly those concerned with the 
nature of the coating, are briefly discussed in the fol- 
lowing section. 


4. NATURE OF THE COATINGS 


Al diffusion coatings, irrespective of the coating tech- 
nique, are characterized by the existence of an alloy-layer; 
in many cases this is comprised of several discrete 
layerstt9), The extent of alloying is dependent on the 
coating process, and the precise nature of this layer is 
discussed later. In addition, part of the coating may be 
of the unalloyed coating metal (e.g. as obtained by con- 
ventional hot-dip processes). The properties of the coat- 
ing will therefore be a complex function controlled by 
the extent and individual properties of the discrete layers 
and their interfaces. 

A cementation coating is a surface alloy-layer of the 
basis and coating metals: it 1s, in most respects, com- 
parable to the alloy-layer of a hot-dipped or annealed 
coating of the same system. Usually cementation coat- 
ings because of their alloy character are harder and less 
ductile than pure metal coatings; they are usually of a 
dull appearance occasioned by the relatively rough sur- 
face produced by the process conditions. Coatings 
produced from the other solid-solid methods also exhibit 


TABLE 4 


Summary of methods for the application of metal coatings 


| Method Comments Eickmess 
(mils) 


Fairly uniform in composition with the exception of some alloys where diffusion at the elec- 
trode surface can result in the deposition of different phases. These coatings are relatively 
thin, and often porous and under stress; they are deposited by discharge from electrolytes. 


Electrodeposition 


eg. Cd, Cu, Cr, Ni, Sn, Zn, Ag, Au, etc. 


Relatively thick coatings, with good coverage, that can be obtained fairly rapidly. The method 


Hot-Dipping 
in. eg. Al, Pb, Sn, Zn, 


is limited to certain metals and the coating tends to be porous or discontinuous if <0.001 


| Good uniformity and close dimensional tolerance can be obtained although coatings tend 


Cementation 


to be slightly thinner than other diffusion methods. There is a limitation on size of article 


| and the coating tends to be slightly harder, but more brittle, because of its alloy nature. e.g. | 


Al, Cr, Zn, etc. 


Cladding 


Metal Spraying 
| operator. e.g. Al, Zn, stainless steel, etc. 


Physical Vapour 


Deposition É 
P compounds can be deposited in this way. 


Chemical Vapour 
Deposition 
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5-10% of substrate 


of oxide. The composition of the coatings is fairly uniform and can be controlled by selec- 
tion of the raw materials. Thick coatings, with good coverage, can be obtained by a skilled 


The method has poor throwing power but can be used for production of decorative electron- 
ic and optical coatings. The main advantage is that all solid elements and thermally stable 


The method has very good throwing power and the coatings are often alloyed to the substrate 
and show good adhesion, e.g. Cr, Mo, Ni, V, etc. 


7 


an alloy-layer but often to a lesser degree and subject 
to the simultaneous influence of other factors. In addi- 
tion to an inner alloy-layer, these coatings also have an 
outer layer of pure coating metal. Hot-dipped coatings 
are similar in this respect but the outer layer is generally 
contaminated with traces of the substrate, an unavoida- 
ble consequence of the mutual alloying action. 
Diffusion can result in the presence of voids and 
porosity(6M. There are, however, instances where inher- 
ent porosity (e.g. the formation of a diffusion coating 
on a sintered-metal substrate) has been decreased by the 
formation of a diffusion coating. Some metals are sen- 
sitive to non-metallic impurities in that they are hardened 
and embrittled by small amounts of carbon, oxygen and 
nitrogen, etc. The formation of coumpounds with these 
elements during diffusion can result in a poorly bonded 
interlayer, a hard or brittle layer, or even a barrier to 
diffusion. As these elements are present in the basis metal 


and sometimes in the process atmosphere, diffusion coat- 


ings frequently exhibit such effects unless special precau- 
tions are observed. Some of these properties can, in fact, 
be used to good purpose. Because of interdiffusion the 
substrate metal can affect the coating markedly; alloy- 
ing elements (as well as non-metallic impurities), in either 
the substrate or coating metal can have a significant in- 
fluence upon the rate of diffusion and structure of the 
coating. In this respect, coatings composed of solid- 
solution series are likely to be more ductile than those 
including brittle intermetallic compounds. 

Generally, there is a measurable dimensional change 
in the coated article, the amount and direction being de- 
pendent on the type and extent of the coating reactions. 
The thickness of alloy-layer can usually be expressed as 
a function of time and temperature although this may 
not always be readily apparent, for example, in hot- 
dipped coatings where crystallites may become detached 
from the alloy-layer and lost in the bath of the molten 
coating metal. 

The presence of internal stresses in coatings is of im- 
portance bacause of their effect on its service life (e.g. 
corrosion resistance, fatigue properties)(6%, Stresses may 
arise from: 

|) Thermal stresses induced at temperatures other 
than the coating temperature and if the coating 
and substrate have different coefficients of ex- 
pansion. 

11) Diffusion stresses leading to volume changes. 

Wi) Anisotropy in the system. 

iv) Surface tension effects caused by the existence of 

preferential substrate sites. 

v) Conversion stresses which may develop when the 
system is taken through a phase-transition tem- 
perature. 

vi) Lattice-misfit stresses caused by epitaxial growth 
(c.f. stresses from anisotropy), occlusions or 
hydrogen absorbed in the lattice. 
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In many instances, the process temperature may be 
sufficiently high to permit stress relief although the ef- 
fect of stress may be reflected in the coating structure 
(e.g. columnar growth). Thick coatings often have a 
strong tendency to crack particularly at corners because 
of the stress introduced during coating formation (c.f. 
oxidation of metals, electrodeposition of metals). Lack 
of space does not allow a further analysis of this sub- 
ject. Supplemented information is provided in Part IH 
of this work. 
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THE DISPERSION EQUATION 
FOR ELECTROSTATIC WAVES IN A BEAM-PLASMA 
SYSTEM IMMERSED IN AN INHOMOGENEOUS 


MAGNETIC FIELD 


C.A.F. Varandas and J.AC. Cabral 


Centro de Electrodinâmica — Complexo Interdisciplinar 


ABSTRACT 


The dispersion equation for electrostatic waves in an 
inhomogeneous beam-plasma system is deduced, in the 
WKB approximation. 


|. INTRODUCTION 


The study of wave propagation in inhomogeneous me- 
dia has great scientific interest because, among other rea- 
sons, inhomogeneities are very commom in natural as well 
as in laboratory plasmas. These inhomogeneities can be 
due to gradients of the plasma density and temperature 
and/or of the confining magnetic field. 

In this work we derive the dispersion equation for elec- 
trostatic waves propagating in a beam-plasma system 
immersed in a practically axial and slowly varying (in z) 
magnetic field. The beam-plasma interaction was one of 
the first subjects to be studied in Plasma Physicst! 2 3), 
Recently the interest in this study results from the possi- 
bility of explaining some effects that occur in fusion 
devices" *º 7) and in ionospheric and magnetospheric 
research'%. In Section IH we deduce the dispersion equa- 
tion for a general stratified medium, defined by its die- 
lectric tensor. In Section HI we obtain the dielectric ten- 
sor and the dispersion equation for the stratified 
beam-plasma system and we analyse a typical dispersion 
diagram. 


2. DISPERSION EQUATION FOR ELECTROSTATIC 
WAVES PROPAGATING IN A STRATIFIED 
MEDIUM FILLING A CYLINDRICAL 
WAVEGUIDE 
The evolution of the electric potential of electrostatic 

waves, propagating in a stratified (in z) mediumf*? fill- 
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ing a cylindrical waveguide, is governed by the Maxwell 
equation 
V eles VO)=0 (1) 
where « is the dielectric tensor of the medium. 
To solve this equation, in cylindrical coordinates, we 
write e under the form”) 


€=6 ElZ) E, AZ) O (2) 
=€, AZ) E(Z) O 
O 0 e, AZ) 


that we will justifify later, and we assume solutions of 
the form 


Q(r, vz D=R(r) Z(z) expl-jmp) exp QGut) (3) 


leading to the following differential equations 


4 
LO dRO qm, 1 1º dO) Lo hã 
R(r) dr” r r Ré(r) dr e 
e, XZ) | d?Z (2) 
elz) Z(z) do 
(5) 
1 del) 1º dZ(7) 
E kZ) dz L(z) dz 


where p is a separation constant that, as we will see later, 
represents the transverse wave number. 
Equation (4) has the general solution 


R(MD)=R, J(pr)+R, Na(pr) (6) 
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where R, and R, are integration constants and J (pr) 
and N (pr) are the m-order Bessel functions of first 
and second kind, respectively. In our case, R, = 0, since 
the electric potential must be finite for r=0, and 


p 
— Fm 7 
ad (7) 
where p,. is the vth zero of J (pr), since the electric 
potential must vanish at the perfectly conducting wall 
of the cylindrical waveguide (r = a). 
Equation (5), if we assume that 


delz) dZ (2) 
dz dz 


dºZ (2) 


<< —+ 8 
dz? (8) 


E, Z) 


can be written under the form 


2 
(O (rm z(=0 (9) 
dz 
where 
E (0) =p? Sat) (10) 
€,, (Z) 


is the dispersion equation for electrostatic waves 
propagating in a stratified medium, filling a cylindrical 
waveguide. 

In the WKB approximation, we can consider zero 
order 


Z(2)=Z, expl-/ k(z) dz] (1) 


and first order 
Z(2) = Z, [k(2)]-12 exp [-/ k(z) dz] 


solutions of equation (9). These two solutions represent 
waves propagating in the positive z-direction, with axi- 


Electron gun Bcrm cavity 


Sd Da pax 


ally varying wave number and either constant amplitude 
(zero order solution) or variable amplitude (first order 
solution). In both cases, the term -["k(2) dz consti- 
tutes, what is sometimes called, the phase memory 
concept, 


3. THE BEAM-PLASMA SYSTEM 


The beam-plasma system (Fig. 1) consists of a cylin- 
drical waveguide, filled with a noble gas at low pressure 
(p = 10") Torr). An electron gun continuously injects a 
beam that creates the plasma by inelastic collisions with 
the neutral gas atoms.-The plasma is confined by a mag- 
netic field that, in this work, we suppose that varies along 
the systems axis. The Maxwell equations lead to the fol- 
lowing radial and axial components of the B-fieldtb 


r cu É d Ba, (2). 
B'(r, z) = 3 E (13) 
2 2 
Bz (r, Z) — B: (z) se A d Ba (z) (14) 


dz? 


These expressions were obtained, for small values of r, 
assuming azymuthal symetry and neglecting the electric 
current density in the plasma. When the magnetic field 
varies slowly on the system's axis [B, (z)], we conclude 
that the radial component B'(r,z) is neglegible and the 
axial component B“(r,z) is practically equal to B. (z). In 
these conditions, the beam-plasma system can be consi- 
dered as a stratificd medium. 


3.1 The dielectric tensor 


The dielectric tensor of the beam-plasma system is 


Fig. 1 Schematic of an experimental beam-plasma system 


(12) defined by('? 
7... DE (SÉ 
did ot e ot (15) 
ES IX ie 
SDS DS St 
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where 


J=En q, (16) 
is the electric current density. In this equation q, 1s the 
charge, n, is the density and v, is the velocity of the 
particles of species a (electrons and ions of the plasma, 
electrons of the beam). 

Assuming that the beam-plasma system can be con- 
sidered as cold (the phase velocity of the wave is much 
greater than the thermal velocities) and as collisionless 
(the wave frequency is much greater than the collision 
frequencies), the movement of each species of charged 
particles is governed byt” 


ip 


dv 2 = 
m = E+v xB 
aa do Se VMB) 


(17) 

The system of equations (15) to (17) can easily be 
solved by a linearization procedure: (i) we consider that 
each dependent variable s (z, t) is separated into two parts: 
an equilibrium part s,(z) that expresses its state in the 
absence of oscillations, and a perturbation part s, (z, t); 
(ii) we neglect the second order terms, assuming that the 
oscillation amplitude is small. 

Considering an uniform neutral plasma and the ions 
as an immobile background of positive charge ensuring 
overall charge neutrality, we have the following fields 


E=E,=-V O, (18) 
B=-B.=B. (zu, (19) 
and we can write 
No, =N,, + O, No, =Np +N di 
Vep = Vip Vip =0 
for the plasma ions and electrons and 
Dep = Do + Dj 
| (21) 
Vob = Vop E Vip 
for the beam electrons. 
The linearization of equation (17) leads to 
o cÉLB Tri | 
m at =-e (E, + B.(Z) Vip X u,) (22) 
for the plasma electrons and to 
(23) 


av av - o 
Be + Va o J==e (E, +B.(g) VeXU) 
for the beam Dl 
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Equations (22) and (23) can be solved, in cylindrical 
coordinates, assuming solutions of the form giving by 
expression (3), leading to the following components of 
the velocities v,, and Vy 


e JwE, —wW.(Z) LM 


vt m— 
P m vaio) vam) 
w.. E + jwE 
VP gs Ce E 9 
di w? -w2. (2) (25) 
É 1! ; 
Via = — a E, (26) 
e JOE ,-u, BE, am 
VI = Edit FEED, = 2 
"om o Gota) (24) 
w. E, + JOE” 
ve Fam Ce l ar Jo 
Ê Qui (2) (28) 
e Ei 
vz pe E | 
Ib m IT? l (29) 
where 
so l d Z(z) ' 
D=ojv, —— ST tY 
Po Za) dz (30) 
The linearized current density is 
J=-e (Mp Vip Elis Va Ft Vob ) (31) 


where we neglect the second order terms n,, v,, and 


Ny Vip: From the beam continuity equation 


V (Mo Vip HM V E + Sato = () (32) 
we can calculate n,, 
l d Z(z) 
Po To do 
js iQ ih Vob (33) 


Substitution of expressions (24) to (33) into equation 
(15) leads to the following results 


o = pe (2) A 
err O Ego O A. — wi. (2) Q —wi.(Z) 58] 
(35) 
= = Se(2) 4 el) copy 
És Li w 2 (2)- ww? Q w” (7)-Q? 


EC 
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2 2 

(2) wns 
= (um 36 
€,; w? Q? ( ) 
€,> a Esr = Es = Cuz = 0 (37) 


where w..(2)/27 is the local electron plasma frequency, 
wp/27 15 the beam plasma frequency and w..(z2)/27 is 
the local electron cyclotron frequency. 

From these results we conclude that the dielectric ten- 
sor of a stratified medium assumes indeed the form given 
by expression (2). 


3.2 The dispersion equation 


Substitution of expressions (34) and (36) into equa- 
tion (10) leads to 


2 op — wplz) 4 
p'|l- 2 : 3 
= We (2) w -w (Z) 


cem [1 DT o o 
e Pio Je 


“that is the dispersion relation for electrostatic waves in 
“a beam-plasma system in a practically axial and slowly 
varying magnetic field. Equation (38) was obtained with 
the first order WKB solution and is only valid when 
both) 

(39) 


= <<] 


É e ER del A 
É od 


2 k'(z) 


| dk (2) 
k? (2) dz” 


and condition (8) are satisfied. For the zero order WKB 
“solution, 


Q=w-k(Z)v, (40) 
“and the dispersion equation is simplified to 
p [1- Vb o Bel?) 
(w-Kk(2)v,,)? cetia, (à Po -w2.(Z) 
, 2 (41) 
+K(g [1- pb Sos 1 =0 
(w- k(Dva) wo? 


When the magnetic field is homogeneous, this equation 
reduces to the well known electrostatic dispersion equa- 
tion for the beam-plasma systemt!2 19, 


3.3 The dispersion diagram 


Fig. 2 shows the dispersion diagram for a set of 
parameters similar to those experimentally determined 
and referred to in our former publicationst!”. For the 
construction of the diagram, equation (42) was numeri- 
cally solved twice, giving either complex-w as a function 
of real-k or complex-k for real-w. In Fig. 2 we can recog- 
nize the electron plasma (a) and cyclotron (b) branches 
of the plasma alone [w,, =O in equation (42)] and the 
two space-charge waves (d and e) and the two cyclotron 
waves (c and f) of the electron beam [w,. = O in equa- 
tion (42)]“2. The complete dispersion diagram presents 
two zones of evanescence (5 and 6), resulting from the 
interaction between the plasma cyclotron branch and the 
beam fast waves, and four instability zones (1 to 4), result- 
ing from the coupling between the plasma modes and 
the beam slow waves. The strongest instability is the ab- 
solute electron cyclotron instability (1). There are another 
absolute cyclotron instability (2), developing at higher 
k-values, and two convective electron plasma instabili- 
ties (3 and 4). 


0 200 


400 600 


FtiMbHz) 


Fig. 2 Dispersion diagram of the beam-plasma system. 
Parameters: p = 1.0cm"!, f. = 350 MHz, f,. = 450 MHz, 


pº[1- o pin, ] 4 fp = 100 MHz, vp = 2.65 X 10ºem s:!. Solid line 
(w-kv ar We wi represents branches of reak-w, real-k solutions: (i) real- 
3 (42) K analysis: Thin solid line is w and dash-dash line is 
2 - 9; (ii) real-w analysis: Dot-dash is k, and dash-dash 
+k Rar ie |- =0 | aço r 
(ww -— ea line is k,. 
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